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■ Abstract 

The existence of stable periodic orbits and chaotic invariant sets of singularly perturbed prob- 
lems of fast-slow type having Bogdanov-Takens bifurcation points in its fast subsystem is proved 
^ I by means of the geometric singular perturbation method and the blow-up method. In particular, 
the blow-up method is effectively used for analyzing the flow near the Bogdanov-Takens type fold 
point in order to show that a slow manifold near the fold point is extended along the Boutroux's 
Q ! tritronquee solution of the first Painleve equation in the blow-up space. 

■ Keywords: fast-slow system; blow-up; singular perturbation; Painleve equation 

1 Introduction 

cn ' 

^ ! Let (xi, • • ■ , jc,j,yi, • ■ • , jm) e R"+'" be the Cartesian coordinates. A system of singularly per- 
\ turbed ordinary differential equations of the form 
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is called a fast-slow system, where the dot ( ' ) denotes the derivative with respect to time t, and 
where e > is a small parameter. Fast-slow systems are characterized by two different time 
scales, fast and slow time. In other words, the dynamics consists of fast motions ({xi, ■ ■ ■ ,x„) di- 
rection in the above system) and slow motions ((yi, • • • , j,„) direction). This structure yields non- 
linear phenomena such as a relaxation oscillation, which is observed in many physical, chemical 
and biological problems. See Grasman [[T3]| . Hoppensteadt and Izhikevich lfT6l and references 
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therein for applications of fast-slow systems. To analyze the fast-slow system, the unperturbed 
system (fast system) of Eq. dl.ll) is defined to be 



ji =0, 



(1.2) 



= 0. 



The set of fixed points of the unperturbed system is called a critical manifold, which is defined 



M = {{xu ■ ■ ■ ,Xn,y\,-- ■ ,ym) e R"^'" I /-(jci , • • • ,Xn,yi,-- ■ ,ym,Q) = 0, ? = 1, • • • (1.3) 

Typically A\ is an m-dimensional manifold. Fenichel IfTTTl proved that if A\ is normally hyper- 
bolic, then the original system (11.11) with sufficiently small e > has a locally invariant manifold 
Ale near A1, and that dynamics on Me is approximately given by the m-dimensional system 



where (xi, - ■ ■ ,Xn,y]_, ■ ■ ■ ,ym) £ R"+'" is restricted to the critical manifold M. The Me is dif- 
feomorphic to M and called the slow manifold. The dynamics of (11.11) approximately consists 
of the fast motion governed by (11.21) and the slow motion governed by (11.41) . His method for 
constructing an approximate flow is called the geometric singular perturbation method. 

However, if the critical manifold M has degenerate points jcq e A1 in the sense that the 
Jacobian matrix df/dx, f = (/i, • • • , fi), x = (xi, ■ ■ ■ , x„) at Xq has eigenvalues on the imag- 
inary axis, then M is not normally hyperbolic near the Xq and Fenichel's theory is no longer 
applicable. The most common case is that df/dx has one zero-eigenvalue at Xq and the critical 
manifold M is folded at the point (fold point). In this case, orbits on the slow manifold Me 
may jump and get away from Me in the vicinity of Xq. As a result, the orbit repeatedly switches 
between fast motions and slow motions, and complex dynamics such as a relaxation oscillation 
can occur. See Mishchenko and Rozov [l25l and Jones IfTSlI for treatments of jump points and 
the existence of relaxation oscillations based on the boundary layer technique and the geometric 
singular perturbation method. 

The blow-up method was developed by Dumortier ^ to investigate local flows near non- 
hyperbolic fixed points and it was applied to singular perturbed problems by Dumortier and 
Roussarie |i7J. The most typical example is the system of the form 



by 



' ji = sgi(xi,-- ■ ,Xn,yi,-- ■ ,ym,0). 



(1.4) 




(1.5) 
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where (x,y) e R^. The critical manifold is a graph ofy = x^ and the origin is the fold point, 
at which the Jacobian matrix of the fast system has a zero-eigenvalue. Indeed, the fast system 
X = -y + x^ undergoes a saddle-node bifurcation as y varies. To analyze this family of vector 
fields, the trivial equation e = is attached as 



with the Jordan block. The blow-up method is used to desingularize such singularities based 
on certain coordinate transformations. The most simple case ^(0,0) is deeply investigated 
by Krupa and Szmolyan et al. [|20l [T2l with the aid of a geometric view point. Straightforward 
extensions to higher dimensional cases are done by Szmolyan and Wechselberger ||33]| for n = 
1 , m = 2 and by Mishchenko and Rozov Il25l for any n and m. Under the assumptions that df/dx 
has only one zero-eigenvalue at a fold point and that the slow dynamics (11.41) has no fixed points 
near the fold point, they show that in the blow-up space, the system is reduced to the Riccati 
equation dx/dy = y - x^ for any n > I and m > 1, and a certain special solution of the Riccati 
equation plays an important role to extend a slow manifold to a neighborhood of the fold 
point, which guides jumping orbits. It is to be noted that the classical work of Mishchenko and 
Rozov [|25l is essentially equivalent to the blow-up method. 

On the other hand, if the dynamics (11.41) has fixed points on (a set of) fold points, for example, 
if g(0, 0) = in Eq. dl.SI) . then more complex phenomena such as canard explosion can occur. 
Such situations are investigated by |l7l|20l[32l|22l|24l| by using the blow-up method. For example, 
for Eq. dl.SI) with g{0, 0) = 0, the original system is reduced to the system x = -y + x^, y = x in 
the blow-up space. If the dimension m of slow direction is larger than 1, there are many types of 
fixed points of (11.41) and thus we need more hard analysis as is done in Il22]| . 

The fast system for Eq. dl.SI) undergoes a saddle-node bifurcation at the fold point. Thus 
we call the fold point the saddle-node type fold point. The cases that fast systems undergo a 
transcritical bifurcation and a pitchfork bifurcation are studied in [|2T]| . It is shown that in the 
blow-up space, systems are reduced to the equations dx/dy = x^ - y^ + A and dx/dy = xy - x^, 
respectively, whose special solutions are used to construct slow manifolds near fold points. 

Despite many works, behavior of flows near fold points at which the Jacobian matrix df/dx 
of the fast system has more than one zero-eigenvalues is not understood well. The purpose of 
this article is to investigate a three dimensional fast-slow system of the form 



X = -y + x^, 
• y = sg(x,y), 
s = 0. 



(1.6) 



Then, the Jacobian matrix at the origin (0, 0, 0) degenerates as 



'0 -1 
^(0,0) 
,0 



(1.7) 



' X = fi{x,y,z,£,6), 
• y = f2(x,y,z,£,6), 
^ z = £g(x,y,z,s,6). 



(1.8) 
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whose fast system has fold points with two zero-eigenvalues, where f\,f2,g are C"' functions, 
e > is a small parameter, and where 5 > is a small parameter which controls the strength 
of the stability of the critical manifold (see the assumption (C5) in Sec. 2). Note that the critical 
manifold 

M(6) = {(X, y, z)eR'\ fy{x, y, z, 0, 6) = f2(x, y, z, 0, 6) = 0} (1.9) 

gives curves on R"* in general. We consider the situation that at a fold point {xo,yo,Zo) e on 
M, the Jacobian matrix d{fi, f2)/d(x,y) has two zero-eigenvalues with the Jordan block, and the 
two dimensional unperturbed system (fast system) undergoes a Bogdanov-Takens bifurcation. 
We call such a fold point the Bogdanov-Takens type fold point. For this system, we will show 
that the first Painleve equation 



dz^ 



,2 =y 



appears in the blow-up space and plays an important role in the analysis of a local flow near 
the Bogdanov-Takens type fold points. This is in contrast with the fact that the Riccati equation 
appears in the case of saddle-node type fold points. It is shown that in the blow-up space, the slow 
manifold is extended along one of the special solutions, the Boutroux's tritronquee solution [[H 
IT9]| . of the first Painleve equation. One of the main results in this article is that a transition map 
of Eq. dl.SI) near the Bogdanov-Takens type fold point is constructed, in which an asymptotic 
expansion and a pole of the Boutroux's tritronquee solution are essentially used. This result 
shows that the distance between a solution of (11.81) near the Bogdanov-Takens type fold point and 
a solution of its unperturbed system is of order 0(s^^^) as £ ^ (see Theorem 1 and Theorem 
3.2), while it is of 0(s^'^) for a saddle-node type fold point (see Mishchenko and Rozov [|25l ). 

It is remarkable that all equations appeared in the blow-up space are related to the Painleve 
theory. For example, the equation dx/dy = y - x^ obtained from the saddle-node type fold point 
is transformed into the Airy equation du/dy = uy by putting x = {du/dy)/u, which gives classical 
solutions of the second Painleve equation. The equation dx/dy = x^ - y^ + A obtained from the 
transcritical type fold point is transformed into the Hermite equation 



d^u _|_ 2 

dy^ dy 



— + 2y— +(A+l)u = 



by putting x + y = -(du/dy) /u, which gives classical solutions of the fourth Painleve equation. 
For other cases listed above, we also see that equations appeared in the blow-up space have the 
Painleve property |l5l[I3; that is, all movable singularities (in the sense of the theory of ODEs 
on the complex plane) are poles, not branch points and essential singularities. This seems to be 
common for a wide class of fast-slow systems. Painleve equations have many good properties 
||5l . For example, poles of solutions of Painleve equations can be transformed into zeros of 
solutions of certain analytic systems by analytic transformations, which allow us to prove that 
the dominant part of the transition map near the Bogdanov-Takens type fold point is given by an 
analytic function describing a position of poles of the first Painleve equation. 



4 



We also investigate global behavior of the system. Under some assumptions, we will prove 
that there exists a stable periodic orbit (relaxation oscillation) if e > is sufficiently small for 
fixed 6, and further that there exists a chaotic invariant set if 5 > is also small in comparison 
with small e. Roughly speaking, 6 controls the strength of the stability of stable branches of the 
critical manifolds. While chaotic attractors on 3-dimensional fast-slow systems are reported by 
Guckenheimer, Wechselberger and Young [[T4]| in the case of n = 1, m = 2, our system is of n = 
2,m = 1. In the situation of [|T4l . the chaotic attractor arises according to the theory of Henon- 
like maps. On the other hand, in our system, the mechanism of the onset of a chaotic invariant 
set is similar to that in Silnikov's works [|28l l29l |30]| . in which the existence of a hyperbolic 
horseshoe is shown for a 3-dimensional system which have a saddle-focus fixed point with a 
homoclinic orbit. See also Wiggins [341 • Indeed, in our situation, the critical manifold M{6) 
plays a similar role to a saddle-focus fixed point in the Silnikov's system. Thus the proof of the 
existence of a relaxation oscillation in our system will be done in usual way: the Poincare return 
map proves to be contractive, while the proof of the existence of chaos is done in a similar way to 
that of the Silnikov's system: as 6 decreases, the Poincare return map becomes non-contractive, 
undergoes a cascade of bifurcations, and horseshoes are created. When one want to prove the 
existence of a stable periodic orbit, it is sufficient to show that the image of the return map is 
exponentially small. However, to prove the existence of a horseshoe, one has to show that the 
image of a rectangle under the return map becomes a horseshoe-shaped (ring-shaped). Thus our 
analysis for constructing the return map involves hard calculations, which can be avoided when 
proving only a periodic orbit. 

Our chaotic invariant set seems to be attracting as that in lfT4l . however, it remains unsolved. 
See Homburg [[T5ll for the proof of the existence of chaotic attractors in the Silnikov's system. 

The results in the present article are used in f3\ to investigate chaotic invariant sets on the 
Kuramoto model, which is one of the most famous models to explain synchronization phenom- 
ena. In [Q, it is shown that the Kuramoto model with appropriate assumptions can be reduced to 
a three dimensional fast-slow system by using the renormalization group method O. 

This paper is organized as follows. In section 2, we give statements of our theorems on the 
existence of a periodic orbit and a chaotic invariant set. An intuitive explanation of the theorems 
is also shown with an example. In section 3, local analysis near the Bogdanov-Takens type fold 
point is given by means of the blow-up method. Section 4 is devoted to global analysis, and 
proofs of main theorems are given. Concluding remarks are included in section 5. 

2 Main results 

To obtain a local result and the existence of relaxation oscillations, the parameter 6 in Eg. (II. 81) 
does not play a role. Thus we consider the system of the form 

' X = fi(x,y,z,£), 

■ y = f2(x,y,z,s), (2.1) 
^ z = eg{x,y,z,e). 
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with C°° functions f\,f2,g'- U x I ^ R, where i7 c R"* is an open domain in and / c R is a 
small interval containing zero. The unperturbed system is given as 

X = Mx,y,z, 0), 

y = f2(.x,y,z,0), (2.2) 
z = 0. 

Since z is a constant, this system is regarded as a family of 2-dimensional systems. The critical 
manifold is the set of fixed point of (12.21) defined to be 

M = {(X, y,z)&U\ Mx, y, z, 0) = f2ix, y, z, 0) = 0}. (2.3) 

The reduced flow on the critical manifold is defined as 

Z = £g(x, y, Z, 0)|(.v,j,,z)eM. (2.4) 

To investigate a Bogdanov-Takens type fold point, we make the following assumptions. 

(Al) The critical manifold M has a smooth component S = S'^^U {L'^} U 5^, where 5^ consists 
of stable focus fixed points, 5,1 consists of saddle fixed points, and where is a fold point. 

(A2) The is a Bogdanov-Takens type fold point; that is, is a Bogdanov-Takens bifurcation 
point of the vector field {fi{x, y, z, 0), fiix, y, z, 0)). In particular. Eg. (12. 21) has a cusp at L^. 

(A3) The reduced flow (12.41) on 5^ is directed toward the fold point and g{L^,0) + 0. 

A few remarks are in order. It is easy to see from (A 1 ) that the Jacobian matrix 5(/i , fi) ld{x, y) 
has two zero eigenvalues at since S,"!" and are saddles and focuses, respectively. Thus there 
exists a coordinate transformation (jc, y, z) ^ {X, Y, Z) defined near such that is placed at 
the origin and Eq. (|2.2l) takes the following normal form 

' X = ai(Z) + a2(Z)Y^ + a3(Z)XY + 0(X\X^Y,XY^, Y^), 
■ Y = bi(Z) + b2(Z)X + 0(X\X^Y,XY^,Y^), (2.5) 
^ Z = 0, 

where ai(0) = bi{0) = so that the origin is a fixed point (for the normal form theory, see Chow, 
Li and Wang [4]). Then the assumption (A2) means that a2(0) 0, a3(0) 0, b2(0) 4^ 0. In 
this case, it is well known that the flow of Eg. (12. 51) has a cusp at the origin (see also Lemma 
3.1). Since Eg. (12. 51) has a cusp at L^, there exists exactly one orbit emerging from L^. The 
assumption (A3) means that if the critical manifold is locally convex downward (resp. convex 
upward ), then g{x,y,z,^) < (resp. g(x,y,z,0) > 0) on 5^^ U {L+}. Thus an orbit of the reduced 
flow on 5^ reaches in finite time. As a result, an orbit of (12.11) may jump in the vicinity of . 
The next theorem describes an asymptotic behavior of such a jumping orbit. 

Theorem 1. Suppose that the system (12.11) satisfies assumptions (Al) to (A3). Consider a 
solution x(0 whose initial point is in the vicinity of 5^^. Then, there exist ?i > such that the 
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distance between x(t), to < t < t\ and the orbit of the unperturbed system emerging from 
is of 0(e4/5) as e ^ 0. 

Note that for a saddle-node type fold points, the distance between x{t) and an orbit emerging 
from a fold point is of To prove the existence of relaxation oscillations, we need global 

assumptions for the system (12.11) . 

(Bl) The critical manifold A\ has two smooth components = U {L^} U and S~ = 
5^7 U {L } U 5,7, where 5* consist of stable focus fixed points, 5* consist of saddle fixed points, 
and where L- are fold points (see Fig|T])- 

(B2) The L- are Bogdanov-Takens type fold points; that is, L- are Bogdanov-Takens bifurcation 
points of the vector field {fi{x,y,z, 0),f2(x,y,z, 0)). In particular, Eq. (|2.2l) has cusps at L-. 

(B3) Eg. (12. 21) has two heteroclinic orbits and a' which connect L^,L~ with points on 5 7, 5^, 
respectively. 

(B4) The reduced flow (12.41) on 5* is directed toward the fold points L- and g{L-,0) 0, re- 
spectively. 

Assumptions (Bl) and (B2) assure that S- are locally expressed as parabolas, and thus they 
are of "J-shaped". Components S ^ and S ~ are allowed to be connected. In this case, 5 ^ U5 is of 
"S-shaped". As was mentioned above, since (12.21) has cusps at L-, there exist two orbits and 
a~ of Eg. (12. 21) emerging from and L . The assumption (B3) means that these orbits are con- 
nected to 5~ and 5^^, respectively. If 5^ U 5" is of "S-shaped", the assumption (B3) is typically 
satisfied because at least near the fold points, the unperturbed system (12.21) has heteroclinic orbits 
connecting each point on 5* to 5*, respectively, due to the basic bifurcation theory. Note that 
the assumption (B3) also determines a positional relationship between 5^ and 5 . For example, 
if 5 ^ is convex downward, S ' should be convex upward. By applying Thm. 1 combined with the 
geometric singular perturbation (boundary layer technigue), we can obtain the following result. 

Theorem 2. Suppose that the system (12.11) satisfies assumptions (Bl) to (B4). Then there 
exists a positive number sq such that Eg. (12. II) has a hyperbolically stable periodic orbit near 
5+ U Qr+ U 57 U Q'" if < e < So- 

To prove the existence of a periodic orbit, the local assumptions are not so important, though 
a positional relationship between components of the critical manifold and the existence of hete- 
roclinic orbits a- are essential. Indeed, similar results for fast-slow systems having saddle-node 
type fold points are obtained by many authors. 

To prove the existence of chaos, we have to control the strength of the stability of 5*. Let us 
consider the system (11.81) with C°° functions f\,f2,g'- t/ x / x /' ^ R, where U and / as above 
and /' c R is a small interval containing zero. The unperturbed system of Eg. (II. 81) is given by 

' X = fi(x,y,z,0,S), 

- y = f2(x,y,z,0,6), (2.6) 
[z = 0. 
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Fig. 1: Critical manifold and the flow of Eq. (l2.1|) with the assumptions (Bl) to (B4). 

The critical manifold M(5) defined by (11.91) is parameterized by 6. At first, we suppose that the 
assumptions (Bl) to (B4) are satisfied uniformly in 6. More exactly, we assume following. 

(CI) There exists Sq such that for every 6 e [0, 6o), the critical manifold Mid) has two smooth 
components 5 +(5) = S^(6) U {L+(5)} US^d) and 5 "(5) = S'id) U {L-(6)} U S;(6). When 6 > 0, 
Sf,(6) consist of stable focus fixed points, 5*(5) consist of saddle fixed points, and L-(6) are fold 
points (see Figll])- Further, the 6 family M{S) is smooth with respect to 5 e [0, 6o). 

(C2) For every 6 e [0,6o), L-(6) are Bogdanov-Takens type fold points; that is, L-{6) are 
Bogdanov-Takens bifurcation points of the vector field (fi(x, y, z, 0, 6), fzix, y, z, 0, 5)). In partic- 
ular, Eg. (12. 61) has cusps at L-(6). 

(C3) For every 6 6 (0, 6q), Eg. (12. 61) has two heteroclinic orbits a^{6) and a'(6) which connect 
L'^(6),L~(6) with points on 5^(5), 5 ^^(5), respectively. 

(C4) For every 6 e [0, 6o), the reduced flow on 5* (5) is directed toward the fold points L-{6) 
and g{L-,0, 6) 4^ 0, respectively. 

In addition to the assumptions above, we make the assumptions for the strength of the stability 
of 5* as follows: 

(C5) For every 5 e [0,6o), eigenvalues of the Jacobian matrix d(fi, f2)/d(x,y) of Eg. (12. 61) at 
(x,y,z) e S^i(6) and at (x,y,z) e S^(6) are expressed by -6 ■ ^t{z,5) + V--Ta»^(z, 5) and -6 • 
IJr{z, 5) + V--Ta»~(z, 6), respectively, where and oj- are real-valued functions satisfying 



The assumption (C5) means that the parameter 6 controls the strength of the stability of stable 



yu*(z,0)>0, a>*(z,0)^0. 



(2.7) 
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focus fixed points on 5* (5). 

Finally, we suppose that the basin of of the unperturbed system can be taken uniformly 
in (5 e (0, 5o): By the assumption (C5), there exist open sets V- D S'^iS) such that real parts of 
eigenvalues of the Jacobian matrix d{fi, f2)/d(x,y) on V- is of order 0(6). In general, the "size" 
of V- depend on 6 and they may tend to zero as 5 ^ 0. To prove Theorem 3 below, we assume 
following. 

(C6) There exist open sets V- d S'^id), which is independent of 6, such that real parts of eigen- 
values of the Jacobian matrix d(fi, f2)/d(x,y) on V- are negative and of order 0(6) as 5 ^ 0. 

This assumption also assures that the attraction basin of of the unperturbed system can 
be taken uniformly in 5 e (0, 6o), see an example below. It seems that (C6) is satisfied in many 
situations because of the existence of heteroclinic orbits a-{6). 

Theorem 3. Suppose that the system (11.81) satisfies assumptions (CI) to (C6). Then, there 
exist a positive number sq and positive valued functions Si(s), Szis) such that if < e < eo and 
6i(s) < 6 < 62is), then Eq.([rHl) has a chaotic invariant set near S+(5) U a^(6) U S~(6) U a-(6), 
where 5i,2(e) ^ as £ ^ 0. More exactly, the Poincare return map 11 along the flow of (11.81) 
near S^i (6) Ua'^ (6) US '^(6) Ua~ (6) is well-defined, and 11 has a hyperbolic horseshoe (an invariant 
Cantor set, on which 11 is topologically conjugate to the full shift on two symbols). 

Theorems 2 and 3 mean that if e > is sufficiently small for a fixed 6, then there exists 
a stable periodic orbit. However, as 6 decreases, the periodic orbit undergoes a succession of 
bifurcations and if 6 gets sufficiently small in comparison with e, then a chaotic invariant set 
appears. In our proof in Sec. 4, 6 will be assumed to be of 0(e(-loge)'''^). We conjecture that 
this chaotic invariant set is attracting, although the proof is not given in this paper. In general, 
given fast-slow systems do not have the parameter 6 explicitly. However, Theorem 3 suggests 
that as s increases for fixed 6, a periodic orbit undergoes bifurcations and a chaotic invariant 
set may appears, see Figl2l Obviously the assumptions (CI) to (C4) include assumptions (Al) 
to (A3) and (Bl) to (B4). In what follows, we consider the system (11.81) with the parameter 6. 
When proving Theorems 1 and 2, 6 is assumed to be constant, and when proving Theorem 3, 6 
is assumed to be of 5 ~ 0(e(-loge)^^^) as e ^ 0. Note that s « e(-loge)'''^ « 1 as e ^ 0. 
Although 5 > is also small, uniformity assumptions on 6 and the fact s « 6 allow us to use the 
perturbation techniques with respect to only on s. 

In the rest of this section, we give an intuitive explanation of the theorems with an example. 
Consider the system 

x = z + 3(y' -y) + 6xi^-y\ 
■ y = -x, (2.8) 
z = esin|^3;j. 

The critical manifold M = A\(6) is given by the curve z = 3(y - y^), x = 0, and the fold points 
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Fig. 2: Typical bifurcation diagram of (11.81) with assumptions (CI) to (C6). 
1 2 

are given by L- = (0, +— — , +— =), see Figl3 

V3 V3 




Fig. 3: Critical manifold of the system (12.81) . 

It is easy to verify that the assumptions (CI), (C2), (C4) and (C5) are satisfied for (12.81) . The 
assumption (C3) of existence of heteroclinic orbits are verified numerically (we do not give a 
proof here). 

The assumption (C6) is also verified by a straightforward calculation. Now we show that 
(C6) implies that the attraction basin of S^(6) of the unperturbed system can be taken uniformly 
in 5 e (0, 5o)- We change the coordinates by an afline transformation (x,y,z) ^ (X, Y,Z) so that 
the point (0, -2/ V3, 2/ V3) is placed at the origin and the linear part of Eq. (|2.8|) is diagonalized. 
Then the unperturbed system of Eg. (12. 81) is rewritten as 
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where the explicit form of the polynomial h, whose degree is greater than one, is too complicated 
to be written here. However, one can verify that h is of the form 



h{X, Y, 6) = ^hiiX, Y, 6) + 6h2{X, Y, 6), (2.10) 

where h\ and h2 are polynomials with respect to X and Y such that all coefficients of hi are 
real. Note that V36 - 5^/2 and 611 correspond to a»^(z, 5) and 5//^(z, 6), respectively, in the 
assumption (C5). 

Now we bring Eg. (12. 91) into the normal form with respect to the first term of the right hand 
side. There exist a neighborhood W of the origin, which is independent of 6, and a coordinate 
transformation (X, Y) i-^ (r, 9) defined on W such that Eg. (12. 91 ) is put in the form 

r = -tr + a^r^ + a5r^ + ■■■ , 
^ e = V36 - 6^2 + O(r^). 

Note that the eguation of the radius r is independent of 6 (see Chow, Li and Wang BUl). In our 
case, as is given by 

-180 + 29^2 

Further, we can prove that a, ~ 0(6), i = 3, 5, • • • as 6 by using the induction together 
with the property that h{X, Y, 0) takes purely imaginary values if (X, Y) e (see Eg. (12. 101) ). See 
Chiba [2] for explicit formulas of normal forms which are convenient for induction. Thus the 
derivative of the right hand side of Eg. (l2.1 II) is calculated as 

± (-^r + a,r' + •■•) = -^(1 + b,r' + hr' + •■■) + 0{6\ (2.13) 

where b^, bs, - ■ ■ are (5-independent constants. It proves that there exists a 5-independent positive 
constant C such that if |r(0)| < C, then r{t) decays as \r{t)\ ~ 0(e'^'^^) for small ^ > 0. The same 
property can be verified for any system with the assumption (C6) by means of the normal form. 

To ascertain the reason why the periodic orbit or the chaotic attractor occur, we take Poincare 
sections Z+,„ 2;„„ S+ and 1+ as in Fig. H 

The section is parallel to the xz plane and located at the right of L^. Take a rectangle R on 
and consider how it behaves when it runs along solutions of Eg. (12. 81) . Since the unperturbed 



system of Eg. (12. 81) has the heteroclinic orbit a+ connecting and S^, the rectangle R also 
approaches to 5^ along and intersects the section Zj^, as is shown in Fig. HI Since the velocity 
s sin(5j/2) in the direction z is positive in the vicinity of 5^ and since 5" consists of stable focus 
fixed points, the intersection area on Ey^ moves upward, rotating around 5^. As a result, the flow 
of R intersects the section E^^, which is parallel to the xy plane, to form a ring-shaped area as is 
shown in Figure IH Further, we can show that the ring-shaped area on Sr moves to 2"^^ along 
solutions of Eg. (12. 81) due to Theorem 1 . The area on 1,'^^^ goes back to the section Z^^^^ in a similar 
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Fig. 4: Poincare sections and a schematic view of the images of the rectangle R under a succes- 
sion of the transition maps. 

manner because Eg. (12. 81) has the symmetry (x,y,z) ^ (-x, -y, -z). Thus the Poincare return 
map n from into itself is well-defined and it turns out that Tl{R) is ring-shaped. 

There are two possibilities of locations of the returned ring-shaped area. If the strength of the 
stability of stable fixed points on 5*, say 6 as in the assumption (C5), is sufficiently large, then 
the radius of the ring-shaped area gets sufficiently small when passing around 5 As a result, the 
returned ring-shaped area is included in the rectangle R as in Figl5](a). It means that the Poincare 
map n is contractive and it has a stable fixed point, which corresponds to a stable periodic orbit 
of Eg. (12. 81) . On the other hand, if the strength 6 is not so large, the radius of the ring-shaped area 
is not so small and it intersects with the rectangle as in Figl5](b). In this case, the Poincare map 
n has a horseshoe. 

(a) (b) /C--^ 

Fig. 5: Positional relationship of the rectangle R with the returned ring-shaped area. 
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3 Local analysis around the fold points 



In this section, we give a local analysis around the fold points Lr by using the blow-up method, 
and calculate a transition map to observe how orbits of Eq. dl.SI) behave near the fold points. To 
prove the existence of chaos, we will give a detailed analysis of the transition map, which does 
not need for the standard proof of the existence of a periodic orbit. The main theorem in this 
section (Thm.3.2) will be made in the end of Sec.3.1. We will calculate only for because 
discussion for is done in the same way. 

3.1 Normal form coordinates 

At first, we transform Eq. dl.SI) into the normal form in the vicinity of L^{5). In what follows, if 
a (formal) power series h centered at the origin begins with n-ih. degree terms {i.e. d'h(0)/dx' = 
(z = 0, 1, ■ • ■ , n - 1) and d"h(0)/dx" ^ 0), we denote the fact as /z ~ Op{n). The notation C>(-) 
is also used to the usual Landau notation. For example if h{x,y,z) ~ 0{x?- ,7^ ,xy,yz,zx) as 
x,y,z^ 0, we simply denote it as h ~ Op(2). 

Lemma 3.1. Suppose (CI), (C2) and (C4). For every 6 e [0,6o), There exists a C°° local 
coordinate transformation (x,y,z) ^ (X, Y,Z) defined near L^(6) such that Eq. dl.SI) is brought 
into the form 

' X = Z - y2 + ciiS)XY + Zhi(X, Y, Z, 6) + Y^hiiX, Y, Z, S) + sh^{X, Y, Z, £, 6), 
< Y = -X + Zh4iX,Y,Z,6) + sh5(X,YZ,s,6), (3.1) 
Z = -£ + sh^iX, Y, Z, e, 6), 

where Ci(5) and /z, (z = 1, • • • , 6) are C°° functions such that Ci(5) > for 5 > and 

/zi, /Z2, /Z4 ~ 0{X, Y, Z), he ~ 0(X, Y, Z, e). (3.2) 

If we assume (C5), then ci(6) ~ 0(6) as 5 ^ 0. 

In these coordinates, L^(5) is placed at the origin and the branch 5 "^(5) of the critical manifold is 
of the form Z = Y^ + C>p(3), X = Op{2). 

Proof of the Lemma. We start by calculating the normal form of the unperturbed system d2.6l) . 
We will use the same notation {x, y, z) as the original coordinates after a succession of coordinate 
transformations for simplicity. Since the Jacobian matrix of (/i,/2) at L^{5) has two zero eigen- 
values due to the assumption (CI), the normal form for the equations of {x,y) is of the form (see 
Chow, Li and Wang [H) 

X = ai{6)z + a2i6)y^ + a3{6)xy + zhi(x,y,z,6) + y^h2ix,y,z,S), 
y = bii6)x + b2(6)z + zh^{x, y, z, S), 

where fli(5),a2(5), (33(5), ^i(5),&2(^) and /zi,/z2,/z4 ~ 0(;c,y,z) are C°° functions. Note that a2(^) 
0, bi(6) for 6 e [0, 6o) because of the assumption (C2). Since we can assume that S ^(6) is 
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locally expressed as z ~ j^, x ~ without loss of generality, by a suitable coordinate transforma- 
tion, we obtain a2{6) = -ai(5) and ^2(^) = 0- Since fixed points on 5 ^((5) are attracting and since 
fixed points on are saddles for 5 > 0, we obtain ai{6)b\{6) < and a^id) > for 5 > 0. 

If we assume (C5), then ajiS) ~ 0(6). We can assume that ai(6) > because we are allowed to 
change the coordinates as x -x, y -y if necessary. Thus, the normal form of Eq. (l2.6l) is 
written as 

X = ai((5)(z + aiid)xy + zhi(x,y, z, 6) + y^h2ix,y, z, 6), 
y = bi(6)x + zh4(x,y,z,S), (3.4) 
^ z = 0, 

with ai(6) > 0, b\{6) < 0. The coordinate transformation which brings Eg. (12. 61) into Eg. (13. 41) 
transforms Eg. (|1.8|) into the system of the form 

' X = ai(5)(z -/) + a3(.S)xy + zhi{x,y,z,S) + y^h2(.x,y,z,6) + £h3{x,y,z,£,6), 
< y = bi{6)x + zh4{x,y,z,S) + £h5{x,y,z,£,S), (3.5) 
^ z = £igi{6) + h(,(x, y, z, £, 6)), 

where h^i^h^^h^ are C°° functions such that h(, ~ 0{x,y,z,£), and where g\{S) := g{V',0,S) is 
a negative constant on account of the assumption (C4). Finally, changing coordinates and time 
scales as 



777 TTTTTTT ' (3-6) 



g,{6)\ ambm 

and modifying the definitions of h'^s (z = 1, • ■ • , 6) appropriately, we obtain Eg. (13. II) . Note that 
since g\{6), ai{5), b\{5) i^0fov6 e [0, 6o), this transformation is a local diffeomorphism for every 
6e[Q,6o). ■ 

Let pi be a small positive number and let 

= {{X, Y,p\) I (X, Y) 6 R2}, Xna = {(X,pI Z) I (X, Z) e R^} (3.7) 

be Poincare sections in the {X, Y, Z) space defined near the origin (see Fig©. The purpose of this 
section is to construct a transition map from St to Recall that there exists an orbit a^{6) 
emerging from L^{6), where L^{6) corresponds to the origin in the (X, 7,2) space. 

Theorem 3.2. Suppose (CI), (C2) and (C4) to (C6). If pi > is sufficiently small, there exists 
eo > such that the foUowings hold for < £ < eo and < 6 < 6^: 

(I) There exists an open set Ue c near the point n 5^(5) such that the transition map 
^/oc • ~^ ^tut along the flow of Eg. (13. II) is well-defined, C" with respect to X and Y, and 
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expressed as 



X 







+ 





(Q + H{X, J/))£^/^ + 0(e log s) , 



(3.8) 



where Q. 3.416 is a negative constant, and Gi, G2, are C°° functions with respect to X, J/, 5. 

The arguments /Y, J/ are defined by 



X = b,{X,Y,Pi,e,5)e-^l' exp[-J(pi, 5)-], 



J/ = D2(X, Y,pus, 6)s-^l' exp[-J(pi, 5)-] 



(3.9) 



where Di^Dj and J are C"" functions with respect to X, Y, 6 such that J > for (5 > 0. Functions 
Di and D2 are not smooth in e, however, they are bounded and nonzero as e ^ and 6^0. 

(II) The point (Gi (pi , 5), p j , 0) is the intersection of a^{6) and S,^,,^. 

(III) The function H satisfies 



dH 

mo) = o, — (^,j/)^o. 



(3.10) 



(IV) If C/g is sufficiently small, for each s g (0, eo) and S g (0, Sq), we can suppose that 

dDi 



dX 

by changing the value of pi if necessary. 

,5^ 



(x,y,pi,£,(5) ^0, 




(3.11) 



Fig. 6: Transition map U^^^ and the heteroclinic orbit a*. 



This theorem means that an orbit of Eg. (II. 81) or Eg. (13. II) running around S^^id) jumps near 
L*i6), goes to the right of L^(6) and the distance of the orbit and the orbit a^(6) is of 0(s'^'^) (see 
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Fig©- In particular, it converges to a'^(6) as e ^ 0. We use the blow-up method to prove this 
theorem. In Sec. 3. 2, we introduce the blow-up coordinates and outline the strategy of the proof 
of Thm.3.2. Analysis of our system in the blow-up coordinates is done after Sec. 3. 3 and the proof 
is completed in Sec. 3. 6. The constant is a pole of the first Painleve equation, as is shown in 
Sec. 3. 3. The function H, which is actually an analytic function, also arises from the first Painleve 
equation. To prove Theorems 1 and 2, it is sufficient to show that X and J/ are exponentially small 
as e ^ 0. However, we need more precise decay rate for proving Theorem 3. For this purpose, 
the factors s'^^^ and s~^^^ will be derived by means of the WKB theory. Eq. (l3.10l) and (13.111) 
are also used to prove Theorem 3. Thus our analysis involves a harder calculation than a usual 
treatment of fold points in fast-slow systems. The assumption (C6) is used to assure that the 
domain Us of the transition map is independent of 5 e (0, ^o)- The assumption (C5) is used to 
show that the argument of exp[- • ■ ] in Eq. (|3.9l ) is of order 0(6). For other parts of the theorem, 
we need only (CI), (C2) and (C4). 



3.2 Blow-up coordinates 

In this subsection, we introduce the blow-up coordinates to "desingularize" the fixed point L^i6) 
having a nilpotent linear part. Regarding e as a dependent variable on t, we rewrite Eq. (l3.1l) as 



f X = Z-Y^ + cii6)XY + ZhiiX, Y, Z, 6) + Y^h2(X, Y, Z, 6) + sh^iX, Y, Z, e, 6), 
Y = -X + Zh^iX, y, Z, 6) + ehiiX, Y, Z, e, 6), 
Z = -£ + sh(,{X, Y, Z, s, S), 
e = 0, 



(3.12) 



with the estimate (13. 2|) . For this system, we define the blow-up transformations Ki,K2 and ^3 to 
be 



{X,Y,Z,e) 
{X,Y,Z,s) 



{r\xi, r\yu r\, r\si). 



iX,Y,Z,£) = (rlx3, rl, /^Z3, rlss). 



(3.13) 
(3.14) 
(3.15) 



respectively, where Ki,K2 and ^3 are defined on half spaces {Z > 0}, {s > 0} and {Y > 0}, re- 
spectively. In what follows, we refer to the coordinates (xi,y\,ri,ei), {x2,y2,Z2, '"2), (-^^3, ri„z3,ei,) 
as Ki,K2,Ki, coordinates, respectively. Transformations Kjj from the Ki coordinates to the Kj 
coordinates are given by 



KU 
K2I 
K32 
[ K23 



t \ t -3/5 -2/5 -4/5 1/5n 

(■^2,j2,Z2,''2) = {XxS^ , yxs^ ' , ' , rxs{ ), 



(xi,3;i,ri,ei) = {X2Z.^''' , y2Z2'\ r2zf, z^""), 

(.X2,y2,Z2,r2) = (.X3£~^'^ 



-1/2 



-2/5 



Z3S- 



,1/4 
2 
-4/5 



-5/4x 



l/5x 

/ \ / -3/2 1/2 -2 -5/2n 

(jC3,r3,Z3,e3) = (X2J2 ' ^2y2 ' ^2^2 ' 3^2 )' 



(3.16) 
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respectively. Our next task is to write out Eq. (l3.12l) in the Ki coordinate. Eqs. (l3.13l) and (13.121) 
are put together to provide 

ii = ri(l -y\ + ci{6)rixiyi + h^{xuyi,ri,6) + y\h<){xi,yi,ri,6) 

3 

+rieihio{xuyurue\,6) + -xie\{l - h-j{xuy\,ri,su5))), 

ji = ri(-xi + rxhn{xuy\,ru6) + r\Eihni.Xuy\,rueu6) + ]:yiSi{\ - h-,{xuyuruSx,8))), 

ri = --r^sx{\ -h-,{xuyuri,ei,6)), 

Si = ^ns^i^ -hj(xuyuri,si,6)), 

(3.17) 

where hj (i = 1, - ■ ■ , 12) are functions such that 

hj{xi,yuri,Si,6) = h(,{rlxurlyur1,rlsi,S), (3.18) 

and hg, ■ ■ ■ ,hi2 are defined in a similar manner through h^, - ■ ■ , /zs, respectively. Thus in these 
functions, xi,yi,Si are always with the factors rj,rj,rj, respectively. This fact will be used in 
later calculations. Note that /z, ~ 0{r^) for i = 7,8,9,11 because of (13.21) . By changing the time 
scale appropriately, we can factor out ri in the right hand side of the above equations: 

' ii = 1 - y^^ + ci(6)rixiyi + h(xuyi,ri,6) + y^ih^ixuyi,^,^) 

3 

+riSihio(xi,yuri,si,6) + -XiSi(\ - hj(xuyi,ri,Si,6)), 

yi = -xi + rihuixuyi,ri,6) + rl£ihu(xuy\,ri,sud) + ^yiei(l - hj(xi,yi,ruSi,6)), 
1 

ri = --riei(l - hj(xi,yuru£u6)). 



(Ki) 



si = ^ei(l -hj(xuyuri,SuS)). 



(3.19) 

Since the time scale transformation does not change the phase portrait of Eq. (l3.17l) . we can use 
Eq. (|3.19l ) to calculate the transition map. 

In a similar manner (i.e. changing the coordinates and dividing by the common factors), we 
obtain the systems of equations written in the K2, Kj, coordinates as 



(K2) 



' X2 = Z2-yl+ r2hi3iX2,y2,Z2,r2,5), 

h = -X2 + rlhu{x2,y2, Z2,r2,6), 

Z2 = -i + rlhi5(X2,y2,Z2,r2,S), 

[r2 = 0, 
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and 



3 2 
;£:3 =-l + Z3+ ci(S)r3X3 + -X2hieixs,r3,Z3,S3,S) + rj^hn(x3,rs,z3,S3,S), 

1 

. r3 = --r3hi(,(x3,r3,Z3,S3,6), 
{1^3) \ I 

Z3 = -S3 + 2z3hi(,{x3,r^,Z3,e3,6) + rle3hi^{x3,r^,Z3,S3,5), 

5 

£3 = -S3hi(,{x3,r3,z3,e3,5), 

respectively, where /zi6(JC3,r3,z3,e3,5) := X3 + rlhi^{x3,r2„z3,S3,5) and hi (i = 13, 
functions satisfying 



(3.21) 



, 19) are 



h„,hi^,hig ~ 0(x3,r3,Z3,S3)- 

Our strategy for understanding the flow of Eq. (l3.1l) near the fold point L*(6) is as follows: In 
Sec. 3. 3, we analyze Eq. (l3.20l) in the Kj coordinates. We will find it to be a perturbed first Painleve 
equation. Since asymptotic behavior of the first Painleve equation is well studied, we can con- 
struct a transition map along the flow of it approximately. In Sec. 3. 4, we analyze Eq. (|3.19l) in 
the Ki coordinates. We will see that in the Ki coordinates, S*(5) has a 2-dimensional attract- 
ing center manifold W^iS) for 5 > (see FigjT]). Since it is attracting, orbits passing nearby 
S'^,(6) approaches W^(6). Thus if we construct the invariant manifold WiS) globally, we can well 
understand asymptotic behavior of orbits passing through nearby 5^,(6). Although usual center 
manifold theory provides the center manifold W'iS) only locally, we will show that there exists 
an orbit y, called the Boutroux's tritronquee solution, of the first Painleve equation in the K2 co- 
ordinates such that if it is transformed into the ^1 coordinates, it is attached on the edge of W(5) 
(see Fig|7]). This means that the orbit 7 of the first Painleve equation guides the manifold W(5) 
and provides a global structure of it. In Sec. 3. 5, we analyze Eq. (l3.21l) in the K3 coordinates. 
We will see that there exists a fixed point whose unstable manifold is 1 -dimensional. Since the 
orbit 7 of the first Painleve equation written in the ^3 coordinates approaches the fixed point, the 
manifold W(5) put on the 7 is also attached on the unstable manifold (see Fig|7]). The unstable 
manifold corresponds to the heteroclinic orbit a'^(6) in the (X, Y,Z) coordinates if it is blown 
down. This means that orbits of Eq. (l3.1l) coming from a region above L^(6) go to the right of 
L*(6) (see Fig© and pass near the heteroclinic orbit q'^((5). Thus the transition map n^^^, is well 
defined. The fixed point in the Kj coordinates corresponds to a pole of the solution 7 in the K2 
coordinates. In this way, the value Q of the pole appears in the transition map (13.81) . 

Combining transition maps constructed on each Kj coordinates and blowing it down to the 
(X, Y, Z) coordinates, we can prove Thm.3.2. 
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Fig. 7: The flow in the (X, Y, Z) coordinates and the blow-up coordinates. The dotted line denotes 
the orbit y of the first Painleve equation. 

3.3 Analysis in the K2 coordinates 

We consider Eg. (13. 201) . Since = e^^^ is a small constant, we are allowed to take the system 



X2 = 22-^2' 
h = -X2, 
Z2 = -1, 



(3.22) 



as the unperturbed system of Eg. (13. 201) . This is eguivalent to the first Painleve eguation : 

d'y2 



dx2 7 

dy2 _ 
—. — - X2, 

\ dZ2 



or 



dzi 



(3.23) 



It is known that there exists a two parameter family of solutions of the first Painleve eguation 
whose asymptotic expansions are given by 

^ N\ ( 1-1/9 lC\ _q/s; _ nZ „ Mr\ . i C2 

^2(22) 

J2(Z2), 



]J-\^2'''-(Y'2'''-'^C2zi'')cOS 

' [ -z5^' + Ciz-'^^cos0 + C2Z2'^'sin0 + O(z2') J' 

(3.24) 
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as Z2 ^ oo and 



^1:2(^2) 

, yiizi) 



^ + |(Z2 - zo) + \{zi - zo)' + 4C3fe - zo)' + 0(iZ2 - zo)') 

(Z2 - zo)^ 5 2 

— ^ + |^(Z2 - Zo)' + ^(Z2 - Zo)' + C3(Z2 - Zo)' + 0((Z2 - Zo)') 
V (Z2 - Zo)^ 10 6 



(3.25) 



as Z2 ^ Zo + 0, where ^ 



4V2 



5/4 



(Z2 



00), and where Ci,C2,C3 and zo are constants 



which depend on an initial value. The value zo is a movable pole of the first Painleve equation 
(see Ince [fTTll . Noonburg [26]. Conte [5]). In particular, there exists a unique solution 7, which 
corresponds to the case Ci = C2 = 0, whose asymptotic expansions as Z2 ^ 00 and as Z2 f2 + 
are of the form 

/ 1 



r 



.^2(Z2) 



(3.26) 



and 



7 



yiizi) 



-12 Q . \ 

— - + -(Z2 - H) + 0((Z2 - nf) 

(Z2 - ii)^ 5 



(3.27) 



respectively, where Q. 3.416. The y is called the Boutroux's tritronquee solution [[TlfT^. 

Let p2 and p3 be small positive numbers and define Poincare sections to be 



|Z2=P2-}, ^"' = {y2=p. 



-2/5, 



(3.28) 



(see FigE]). By Eqs. (l3.26l[X27]) . the intersections P2 = yPiS^"', Q2 = yCM.^ of y and the sections 
are given by 

P2 = (p., Py, Pz) = [-{2l3fl^pfl' + 0(pf ), pfl\ a + V6pJ^' + 0(p3)) , (3.29) 
Q2 = (q., qy, qz) = (-pf 1'^ + 0(pf'\ -p'^" + 0{pf), p-^'") , (3.30) 

respectively. 

Proposition 3.3. If p2 and p3 are sufficiently small positive numbers, there exists an open set 
U2 c such that the transition map 112'^ : U2 ^long the flow of Eq. (l3.20l) is well-defined 

and expressed as 



n 



loc 



' X2 ' 








yi 




Py 

Pz 


+ 


K ^2 ; 




lo J 


V 





ii2{x.2 - qx,y2 - qy,P2, r2,P3, 5) 
r2 



(3.31) 



where Hi{x,y,p2,r,pT„6) and H2{x,y,p2,r,p2,6) are C°° functions with respect to x,y,r and 6 
satisfying the equalities //i(0, 0,p2, 0,p3, 5) = H2(0,0,p2,0,p3,6) = for any small p2,P3 > 
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Fig. 8: The solution y of the first Painleve equation and the Poincare sections, 
and 6 e [0, ^o). 

Proof. This is an immediate consequence of the differentiability of solutions with respect to 
initial values X2,y2 and parameters r2, 6. Note that at this time, we did not prove differentiability 
at p3 = 0, which will be proved in the next Lemma. ■ 

Since Hi and H2 are C°° with respect to r and 6, we put them in the form 

Hi{x,y,p2,r,p3,6) = Hi(x,y,p2,P3) + 0(r), z = 1,2, (3.32) 

where we use the fact that when r2 = 0, the system (13.201) is independent of 6. Then, the value 
limpj^o [pz + H2(x - qx, y - qy,P2,p^)) gives a pole zo of a solution of Eq. (|3.23l) through an initial 
point {x,y,p'^^'^)\ that is, X2(z2),y2fe) ^ 00 as Z2 ^ Zq. Prop. 3. 3 implies that Hi are in x and 
y when p^ > 0. Now we show that Hi can be expanded in p^^^ and they are C°° even if p3 = 0. 
This means that a position of a pole is also smooth with respect to initial values. In the proof, the 
Painleve property will play a crucial role. Part (ii) of the next Lemma is used to prove Thm.3.2 
(III). 

Lemma 3.4. (i) The functions Hy and H2 are analytic with respect to {x,y) 6 U2, P2^ > and 
P3^^ > 0, though they are singular at p]^^ = 0. 

(ii) ^2(0,0,p2,0) = 0, ^H2{x,y,p2,Q)^Q. 

ox 

Proof. Let X2 = X2{z2, Pi, xq, yo) and 3^2 = 3^2(^2; P2, xq, yo) be a solution of the system (13.231) with 
the initial condition 

X2{p~2^'^;p2,XQ,yo) = xo, y2(P2'^'^;p2,xo,yo) = yo. 

Suppose that y2iz) = p^^^^ for some z = zixQ,yQ,p2,P3). When p3 > 0, the statement (i) imme- 
diately follows from the fundamental theorem of ODEs: Since the right hand side of the system 
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(13.231) is analytic, any solution is analytic in time zi, initial time p^^'^ and initial values {xQ^yo). 
Applying the implicit function theorem to the equality 



-2/5 



one can verify that 



Z(Xo,yo,P2,p3) = Pz + HiiXo - qx,yo - qy,p2,P3) 



is analytic in xo,yQ,P2^ > and pl^'^ > 0. Thus 



X2(z(xo,yo,p2,P3y,P2,XQ,yo) = Px+Hi{xo - qx,yo - qy,P2,P3) 



(3.33) 
(3.34) 

(3.35) 



is also analytic in the same region. Since z ^ oo as p2 ^ 0, Hi and H2 are singular at p^^ = 0. 

When p3 = 0, z(xo,yo,p2, 0) gives a pole and ^2 = = 00 at Z2 = z(xo, jo,P2, 0). Thus we 
should change the coordinates so that a pole becomes a regular point. For (13.231) . change the 
dependent variables (^2, J2) and the independent variable Z2 to (^, rj) and r by the relation 



T]^ 2 2 

y2 = — J, 



(3.36) 



and Z2 = KT, respectively, where k := (-6)^^^ < 0. Then, (|3.23l) is brought into the analytic 
system 



1 



Jr 8 



8 



1 



,3 5 4 ^ 



(3.37) 



Since any pole of 3^2(22) is second order [fTTll . a pole of );2 is transformed into a zero of 77(7) of 
first order. Let 77 = //(t; 5, 770, ^0) and ^ = ^(t; 5, 770, ^0) be a solution of the system satisfying the 
initial condition 



77(5; s, 770, ^0) = .0, ^(*; .0, ^0) = ^0, 

where (770, ^0) and the initial time s correspond to (xcjo) and p^"^^^, respectively, by the transfor- 
mation (13.361) . Suppose that 

77(t(5, 770, ^0,P3);*, .0,^0) = i-K^y^pl'^ 

for some r = f(s, 770, ^cPs), which corresponds to a value of z(xo,yo,p2,P3) by the relation z = kt 
so that 3^2(2) = P3^^^ (note that when 3^2 = P3^^^ then 77 = {-K^y^p^J^). Since 



dr] 
57 



;,=(-^3)l/2^>/5 



= 1 + Oipf), 
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the implicit function theorem proves that f is analytic in s, /7o>^o and small p^^^ > 0. Since the 
transformation (770, ^0) ^ i^o, yo) defined through (13.361) is analytic when yo 4^ 0, it turns out that 
z(xo,yo,P2,P3) is analytic in (jco,yo) e U2,P2^ > and pl'^ > 0. Now Eqs. (l3.34l 13.351) prove the 
part (i) of Lemma. 

To prove (ii), let us calculate the asymptotic expansion of f(s, 770, ^0, 0), at which rj = 0. We 
rewrite (13.371) as 

^ dT 1 



du~ l + lr]^ + |77' - \7i^^' 

d^ |t^?7 + W - (^^ - i) - 1^'^ + 

A general solution of this system is obtained in a power series of r\ as 



(3.38) 



^15 



1 



^17 



T = Tx+r]- ^if - —rf + ^n' + Oirf), 
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12 



14 



5ti 



16 



24 



(3.39) 



where T\ and ^1 are constants to be determined from an initial condition. By using the initial 
condition (r, 77, ^) = {s, r]o, ^0), ti is determined as 



1 



^0 7 



(3.40) 



When ?7 = 0, T = Ti . This means that the above t\ gives the expansion of f{s, 770, 0). Then we 
obtain 



^^2, dz 
-r—{xo-qx,yo-qy,P2,0) = —{XQ,yo,p2,0) 
OXq OXq 

df 

= K—{s,r]o,^Q,Q) 
0x0 

df drjo df d^Q 
= K- — h a:- 



5770 dxo d^o dxo 



3 ' 



which is not zero for small 770 (thus for large yo). The equality ^2(0, 0,p2, 0) = is obvious from 
the definition. ■ 



Remark. Since //, is analytic in > 0, it is expanded as 



Hi(x,y,p2,P3) = Hi(x,y,p2) + 0(p\'^), 



(3.41) 
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for / = 1,2. Indeed, one can verify that 

3 V6 

Hi{x,y,p2,P3) = Hi(x,y,p2,0) + V6pf + —(Ht(x,y,p2,0) + p,)p, + 3pf + 0(pf ) 

by using the expansion (13.251) . Further, Hi are expanded in a Laurent series of p^^. In particular. 
Eg. (13. 401) show that the expansions are of the form 

Hi(x,y,p2) = ^i(x,y) + p-''' Fi(x,y,p'^"\ (3.42) 

because s = Pj'^^^ /k, where Fi,F2 are analytic functions. The proof of the above lemma is based 
on the fact that a pole of (13.231) can be transformed into a zero of the analytic system by the 
analytic transformation. This property is common to Painleve equations, and the transformation 
(13.361) is used to prove that (13.231) has the Painleve property llSlfTTll. 

3.4 Analysis in the Ki coordinates 

We turn to Eg. (13. 191 ). It is easy to verify that Eg. (13 . 1 91) has fixed points (xuyi,ri,si) = (0, ±1,0, 0) 
By virtue of the implicit function theorem, we can show that there exist two sets of fixed points 
which form two curves emerging from (0, +1,0, 0), and they correspond to 5^(5) and re- 
spectively (see FigjV]). On the fixed points, the Jacobian matrix of the right hand side of Eg. (13. 191) 
has eigenvalues given by 

0, 0, ^ {cMnyi + 0(rl) ± ^Sy, - 4ci{6)r,xi + ^(r^)) . (3.43) 

In particular, the eigenvalues become 0, 0, + V2/ at the fixed point Qi = (0, -1, 0, 0), but at fixed 
points in S'^(6)\Qi, they have two eigenvalues whose real parts are negative if ri is small and 
S > 0. Eigenvectors associated with the two zero eigenvalues at points on 5^^(5)\2i converge to 
those at Qi, which are given by (0,0, 1,0) and (-1,0,0,2), as 0. The vector (0,0, 1,0) is 
tangent to S'^{6). Thus (-1,0, 0, 2) is a nontrivial center direction. 

Lemma 3.5. If 5 > 0, there exists an attracting 2-dimensional center manifold W^id) which 
includes 5^(5) and the orbit y of the first Painleve eguation written in the Ki coordinates (see 
Fig©. 

Proof. Let B(a) be the open ball of radius a centered at Qi. Since at points in 5'^^(5)\5((3) the 
Jacobian matrix has two zero eigenvalues and the other two eigenvalues with negative real parts, 
there exists an attracting 2-dimensional center manifold 1^(5, a) emerging from S^{6)\B{a) for 
any small a > 0. Let y be the solution of the first Painleve eguation described in the previous 
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subsection. Its asymptotic expansion (13.261) is written in the coordinates as 



yi 



-1 + 0{zf'^) 


-5/4 



(as Z2 °°), 



(3.44) 



by the coordinate change K2i (13.161) . The curve (13.441) approaches the point as Z2 ^ oo and 
its tangent vector converges to the eigenvector (-1,0,0,2) at Q\ as zi oo. Thus W\6) : = 

lim W^{6, a) U y forms an invariant manifold. ■ 

«— >o 

Note that y is included in the subspace {r\ = 0}. This lemma means that the orbit y guides 
global behavior of the center manifold W((5). 

Let pi,p2 > be the small constants referred to in Thm.3.2 and Prop. 3. 3, respectively. Take 
two Poincare sections E"' and S™" defined to be 



S'l" = {(jci,yi,ri,ei)|ri = pi,\xi\ < pi, \yi + 1| <pi, < ej <p2}, 
{(xuyuri,si)\0 < n <pi,|xi| <pi, \yi + 1| <pi, si = pi). 



(3.45) 



respectively. Note that E"' is included in the section St (see Eg. (13. 71) ) if written in the (X, F, Z) 
coordinates and Sj"' in the section (see Eq. (l3.28l) ) if written in the K2 coordinates. 




Fig. 9: Poincare sections to define the transition map 11^. 



Proposition 3.6. Suppose (CI), (C2) and (C4) to (C6). 

(I) If pi and p2 are sufficiently small, the transition map n'°^ : 2"' 1,""' along the flow of 
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Eq. (|3.19l) is well-defined for every 6 e (0, 60) and expressed as 



loc 



Xi 




( I 1/5 -1/5 cnN 






yi 




t 1/5 -1/5 c\ 


+ 




pi 




1/5 -1/5 
Pl^l P2 





.^1 , 




P2 




I J 



(3.46) 



where ^1 and are C°° functions such that the graph of x\ = ipi{ri,si,6) and yi = (p2{ri,si,6) 
gives the center manifold W^id). The second term denotes the deviation from W^iS), and and 
Yi are defined to be 



/P2\ r ^1 
Xi = Di(xuyuPi,£i,P2,5)\ — \ exp[-(i(pi,ei,p2,5)— J, 

Fi = D2(xi,3;i,pi,ei,p2,(5)l — I exp[-J(pi, ei,p2, 5)— J, 



(3.47) 



where £>i, £>2 and J are C°° functions with respect to xi,yi,pi and 6. Although D^Dj and d are 
not C"' in £1 and p2, they are bounded and nonzero as £1 ^0 and 6^0. Further, they admit the 
expansions of the form 



Di(xuyuPuSup2,S) = DiixuyuPi,SuS) + 0((ei/p2)'^^), 
d{pusup2, 5) = d{pi,5) + 0((ei/p2)^^^), 



(3.48) 
(3.49) 



for i = 1,2. 

(II) The first term in the right hand side of Eg. (13. 461) is on the intersection of Sj"' and the center 
manifold W^iS). In particular, as s\ 0, n'"'(jci, ji,pi, eO converges to the intersection point of 
2';"' and 7. 

(III) If the initial point {x^ , , pi , eO is sufficiently close to W\6), 



dx 



-{xuy\,Pu£\,5) + 



(3.50) 



except for a countable set of values of S\ . 

Remark. To prove the existence of a periodic orbit, it is sufficient to show that X\ and Y\ are 
exponentially small as s\ 0. However, to prove the existence of chaos, we need more precise 
estimate as the factors {p2l^\9''^ and (p2/ei)^^^. Eg. (13. 501) is used to prove Eg. (l3.11l) . 

Proof. At first, we divide the right hand side of Eg. (13. 191) by 1 - /z? and change the time scale 
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accordingly. Note that this does not change the phase portrait. Then we obtain 

3 

xi = 1-)^]+ ci(6)rixiyi + -x^si + h+ y\h<) + riSihiQ 

+(1 -y\+ ci{6)rixiyi + h+ y\h<) + riei/zio)/z2i, 

ji = -xi + ^yisi + rihn + rlsihu + {-xi + n/zn + r\eihn)h2i, (3.51) 
1 

ri = --^riEu 

■ _ ^ 2 

£1 - ^Si, 

where /z2i = ^^d arguments of functions are omitted. Equations for ri and Si are solved 

as 

ri(0 = ri(0)^ j , .,(0 = ^-3^, (3.52) 

respectively. Let T be a transition time from 2"' to Sj'"^. Since S\{T) = p2, T is given by 

T'll-m (3.53) 



5ei(0) \ p2 

To estimate xi(T) and yi(T), let us introduce the new time variable r by 



4-5gi(0)f j 



1/5 

(3.54) 



Then, ri(f) = ri(0)T, £i(0 = ei(0)T ^. Note that when t = 0, t = 1 and when t = T, one has 

Claim 1. Any solutions (xi,yi) of (13.511) are of the form xi = r~^u\(T), yi = t~^U2{t), where ui 
and U2 are C"' with respect to t. 

Proof. Changing the time t to r, the system (13.511) is rewritten as 

1 3 
-Tei(0)T-4-^ = l-yl+ c,(6)rmrxm + T^iei(0)T-^ + h,+y]h, + ri(0)ei(0)T-^/zio 
4 4 

+{l-yt + ci{6)ri(0)TXiyi + h+y% + ri(0)ei(0)T %w)h2u 

-7ei(0)T-4^ = -^1 + i^iCO)^"')'! + n(0)T/z„ + niQfsmr-'hn 
4 ar 2 ^ , 

+(-Xi + rMThn + r,{0femT~^hn)h2i. 

Putting x\ = T~^u\, y\ = t~^U2 yields 

-^ei(O)^ = (r^ - T^ul + ci((5)ri(0)r^MiM2 + T^h + r^w^/zg + ri(0)ei(0)T^/iio)(l + /zii), 
t (3.55) 



4 ' dr 

11 



Recall that /z? is defined through (13.181) . and thus 

h,{x,,yurusu5) = h(rm'uun{0fu2,rM'r\rM'sm,d), (3.56) 

which implies that is C" with respect to mi, M2, ri(0), ei(0), 5 and r. Functions h^, ■ ■ ■ ,h\2 
and /i2i have the same property. Hence the right hand side of Eg. (13. 551) is C°° with respect to 
Ml, M2, ''i(O), 6 and t, which proves that solutions U\{t) and M2(t) are C°° with respect to ri(0), 5 
and r. ■ 
Next thing to do is to derive the center manifold and how x\{t) and y\{t) approach to it. 
The local center manifold W^id) is given as a graph of functions x\ = (pi(ri,si,6), yi = 
<P2(f\,Si,S). By using the standard center manifold theory, we can calculate ipi and ^2 as 

(piiruSi,S) = -^£1 + 0{rl,ri£i,sl), (p2iruSi,S) = -1 + 0{rl,riSus])- (3.57) 

To see the behavior of solutions xi and yi near the center manifold W^id), we put xi and yi in the 
form 

^i(T) = ^i(ri(T), ei(T), 6) + t-\i(t), jKt) = ^2(^1(7), ei(T), 5) + t-2v2(t). (3.58) 

Since t^xi(t) and t^Ji(t) are C°° in t for every solutions xi and ji, so are solutions T^(pi(ri(T), si(t), 6) 
and T^(p2iri(T), siir), S) on the center manifold multiplied by and r^, respectively. This implies 
that vi(t) and V2(t) are also C°° in r. Substituting Eg. (13. 581) into (13.511) and expanding it in vi, V2 
and ei(0), we obtain the system of the form 

Si-p- = -8t^V2 +4ciriT^Vi + r\T^h22{ry,T,5)vi + r\T^ h2^{ri,T,5)V2 + ^i(vi, V2, Sj, 5, t), 

JV2 ^^-^^^ 

ei— = 4t\i + r|T^/?24(ri, T, 5)vi + rlr' h25(ri,T, 6)v2 + ^2(vi, V2, ri, Si, 5, t), 

where ^1,^2 ~ 0(vl,ViV2,vl,£i) denote higher order terms, h22, • • • ,/z25 are C°° functions, and 
where ri(0),£i(0) and ci(6) are denoted by ri,ei and Ci, respectively. This is a singular perturbed 
problem with respect to si . 

Claim 2. Any nonzero solutions of this system are expressed as 

vi = Di(T,ri,ei,5;vio,V2o)exp , V2 = Z)2(t, rj, £1, 5; Vio, V20) exp , 

(3.60) 

where Vio = Vi(l) and V20 = V2(l) are initial values, and where and d* are in 

T, ri,vio,V2o and 6. Although D* and D* are not in £1, they are bounded and nonzero as 
si ^ 0,6 ^ 0. If vio, V20, n and t are sufficiently small, 

dD\ 

- — (T,ri,ei,5; vio,V2o) 0, (3.61) 

OVw 
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except for a countable set of values of Si . 

Proof. At first, we consider the linearized system of (13.591) as 



dvi 
dr 



s\ — = -8r^V2 + AariT^y + r\T^h22{ruT,6)vY + r^r'^ h23(ri,T,6)v2 + <9(ei), 

dV2 ^^-^^^ 
dr 

which yields the equation of Vi as 

- sMcir.r' + 4ryh26 + 0(e,))^ + (32t^ - 4r?r^/i27 + 0(si))vi = 0, (3.63) 
dT'^ dr 

where h2(,{r\,T, 6) and /?27(ri, t, 5) are C°° functions. According to the WKB theory, we construct 
a solution of this equation in the form 



Vi(t) = exp[— ^ ei5„(T)]. 



"^1 «=() 

Substituting this into Eq. (l3.63l) . we obtain the equation of 5o(t) 

- (4ciriT^ + Ar\T'h2/4^ + 32r' - Ar\T^h2, = 0, 
dT I dr 

which is solved as Sq = Sq{t) = V{t) + iW{T), where 



J^ilcyris^ + 2r\s^h2e)ds, W{t) = J^(. 



V(t)= I (Icons' +2ris'h2e)ds, W(t) = ( (2ciri5^ + 2r>'/z26) a/: . 3 ^^ -Ids, 



Ss^ - r\s''h2i 



are real- valued functions for small ri. If ri > is sufficiently small and if ci(5) > 0, < t < 1, 
then V{t) < 0. For these 5 q (r) and 5 q (t), S *(r), 5 2 (t), • • • are uniquely determined by induction, 
respectively. Thus a general solution vi(t) is of the form 

vi(t) = r exp[y(T)/ei] exp[/W^(T)/£i] exp[5| + + + • • • ] 

-l-^"exp[V(T)/ei]exp[-zW(T)/ei]exp[57 + S\S2 +■•■], 

where ,kr eC are arbitrary constants. Put 

= exp[iW(T)/si] exp[5[ + siS^ + ■■■], D\_ = exp[-?W(T)/ei] exp[5r + S1S2 +■■■]■ 

Then, vi is rewritten as 

vi(t) = r exp[\/(T)M]Z)*^ +r exp[V(T)/ei]Z)*_, 
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where and D\_ are in Vio, V20, r, and 6. They are not in £1 because of the factor 
1 /ei, however, they are bounded and nonzero as s\ ^ 0. In a similar manner, it turns out that V2 
is expressed as 

V2(t) = r exp[y(T)/ei]D*^ + k~ exp[y(T)/£i]D*_, 

where D*^^ and are C°° in vio, V20, t, r\,5, and are bounded and nonzero as £1 0. Therefore, 
the fundamental matrix of the linear system (13.621) is given as 

F{T) = I j exp[y(T)/ei]. (3.64) 

Now we come back to the nonlinear system (13.591) . We rewrite it in the abstract form as 

dv 

£1— =A(T)v + g(v,T), 
dr 

where v = (vi, V2), g = igi,g2), and A(t) is a matrix defining the linear part of the system. To 
estimate the nonlinear terms, the variation-of-constants formula is applied. Put v = F(t)c(t) with 
c(t) = (ci(t), C2(t)) e C^. Then, c(t) satisfies the equation 

^ = -F(Tr'g(F(T)c,T). (3.65) 
dr £1 

Let c = c(t, £1) be a solution of this equation. Since F(t) ~ 0(e^^^^^'^') tends to zero exponentially 
as ei ^ and since g is nonlinear, the time-dependent vector field defined by the right hand side 
of (13.651) tends to zero as ei 0. Since solutions c(t, Si) are continuous with respect to the 
parameter £1, it turns out that c(t, Si) tends to a constant as ei 0, which is not zero except 
for the trivial solution c(t, ei) = 0. This proves Eq. (l3.60l) with the desired properties by putting 
d* = -V{t) and D* = D^ci + D*_C2 {i=\,2). Note that since the right hand side of (l339l) is not 
zero at 5 = 0, ^ 0,D* ^ as 5 ^ 0. 

When ri = vio = Vjo = 0, the derivatives dvi/dvio, (i = 1, 2) with respect to the initial value 
Vio satisfy the initial value problem 



£y — ^iT,0,£,,6;0;0) = -ST'-^iT,0,£uS;0;0), -^(l,0,sy,6;0;0) = 1, 

dr dvio dvio ovw 

si 3- ^(r, 0, £i , 6; 0; 0) = 4t' -^(r, 0, £1 , 5; 0; 0), -^{1,0, £u 6; 0; 0) = 0. 

dr dviQ dvio dvw 



(3.66) 



This is exactly solved as 



5vi f4V2 , 

-^(r, 0, £uS; 0; 0) = cos — (r^ - 1) 
dvio \ 5£i 



(3.67) 
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In particular, 

1^(0, 0, Si, 6; 0; 0) = cos [i^] (3.68) 

is not zero except for a countable set of values of si. This and the continuity of solutions of ODE 
prove Eq. (l3.61l) . ■ 
Let us proceed the proof of Prop. 3. 6. For Vi(t) and V2(t) in (13.601) . Xi(t) and yiir) are given 
as (13.581) . Since r = (ei(0)/p2)'''^ when t = T,we obtain 

+ Z)i((ei(0)/p2) ' , ri(0), £,(0), 5; Vio, V20) exp — , 

\ei(0)/ L ei(0) J 

3'i(r) = V2(rl(0)(el(0)/p2)^/^P2,5) 



+ 



— — D2((ei(0)/p2)''', ri(0), ei(0), d\ Vio, V20) exp — 

ei(0)/ L £j(0) 



Put 

D;((el(0)/p2)^^^ ri(0), ei(0), 6; vio, V20) = A(^i(0),);i(0), ri(0), ei(0),p2, 5) 

for z = 1,2. Since D* is C°° in vio, V20, niO) and 6, Di is also C~ in xi(0),yi(0), ri(0) and 5. Since 
D* is C°° in t, is bounded and nonzero as ei(0) 0. Finally, let us calculate 



d*((£m/P2r", rm, S) = (IciiSymr' + IniOrr'hjMO), T, S))dT. 

'(£l(0)/p2)'/5 



Due to the mean value theorem, there exists a number r* > such that 



d*((em/P2y",rm,S) = ^ci(6)rM 



1 



em 

P2 



+ /Z26(ri(0),T*,5)^^ 



By the assumption (C5), an orbit of (13.511) near the center manifold W^iS) approaches to 1^(5) 
with the rate 0{e~^^*'). By the assumption (C6), such an attraction region (basin) of W^id) exists 
uniformly in 5 > at least near the branch 5^^(5). Thus /i26('"i(0)^ t*^ ^) is of order 0(6) as well as 
Ci(S) if P2 > is sufficiently small. Therefore, there exists a function d, which is C°° with respect 
to ri(0) and 6, such that 

d*((sm/P2)''\n(0XS) = d(rm,em,P2,S) ■ 6. 

Since p+(z, 0) ^ 0, J(ri(0), ei(0),p2, 0) ^ 0. Since D* and d* are in r = (el/p2)^^^ they admit 
the expansions (I3.48[ 13.491) . This proves (I) of Prop. 3. 6. Proposition 3.6 (U) is clear from the 
definition of (pi,(p2, and (III) follows from Eg. (13. 611) . ■ 
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3.5 Analysis in the Ki, coordinates 

We come to the system (13.211) . This system has the fixed point (X3, r^, Z3, £3) = (- V2/3, 0, 0, 0) 
(see Fig|7]). To analyze the system, we divide the right hand side of Eg. (13. 211) by -hi^(x3, r3, Z3, £3, 6) 
and change the time scale accordingly. Note that this does not change the phase portrait. At first, 
note that the equality 



1 



hi(,(X3,r3,Z3,S3,5) 



+ h3i(x3+ yf2/3,r3,z3,S3,5) 



(3.69) 



holds, where /?3i ~ Op{2) is a C°° function. Using Eg. (13. 691 ) and introducing the new coordinate 
by X3 + V2/3 = X3, we eventually obtain 



X3 = -3x3 + 

1 

Z3 = -2z3 - 



S3 - "2^3, 



-Z3-Ci(6)r3 + h32(X3,r3,Z3,S3,6), 



-£3 + S3h33(X3,Z3,£3,S) + S3r3h34iX3, r3, Z3, £3, S), 



(3.70) 



where h32 ~ Op(2) and /i33,/z34 ~ Op{\) are C°° functions. Note that /?33 is independent of r3. 
This system has a fixed point at the origin, and eigenvalues of the Jacobian matrix at the origin 
of the right hand side of Eg. (13.701) are given by -3, 1 /2, -2, -5 /2. In particular, the eigenvector 
associated with the positive eigenvalue 1/2 is given by (-2ci((5)/7, 1,0,0) and the origin has a 
1 -dimensional unstable manifold which is tangent to the eigenvector. The asymptotic expansion 
(|3.27l) of the solution y of the first Painleve eguation is rewritten in the present coordinates as 



{X3,r3,Z3,S3) = {0((Z2 ' ^T), 0, 0{{Z2 ' ^f), 0((Z2 ' n)^)). 



(3.71) 



which converges to the origin as Z2 ^ (see FigfTOl). 

Let pi and p3 be the small constants introduced in Sec. 3.1 and Sec. 3. 3, respectively. Define 



Poincare sections X'" and X""' to be 



^3 = {(X3,r3,Z3,S3)\\X3\ <Pl, < r3 <Pi, IZ3I <pi, £3 =P3}, (3.72) 
23"' = {iX3,r3,Z3,S3)\\X3\ < Pl, = pu fcl <Pi, < £3 < P3}, (3.73) 



respectively (see FigfTOl). Note that S"' is included in the section SS"' (see Eg. (13. 281) ) if written in 



the K2 coordinates and S""' in the section (see Eg. (13. 71) ) if written in the (X, Y, Z) coordinates. 
Proposition 3.7. (I) If pi and p3 are sufficiently small, the transition map 113'"' : S3"' 
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Fig. 10: Poincare sections to define the transition map 11; 



loc 



along the flow of Eq. (l3.70l) is well-defined and expressed as 



T-r/oc 
^^3 



\P3j 



Pi 



(3.74) 

where /3i and /S^ are C°° in their arguments, and are C°° with respect to x^, Z3,P3 and 6 with 
the property that /32 and /34 are bounded as ri, 0. 

(II) As r3 ^ 0, n3''^(^3, r3,z3,p3) converges to the intersection point (/3(pi,5),pi,0,0) of S""' 
and the unstable manifold of the origin. 

Before proving Prop.3.7, we need to derive the normal form of Eg. (13. 701) . 

Lemma 3.8. In the vicinity of the origin, there exists a C°° coordinate transformation 



r3 

Z3 
^3) 



= 0(X3,r3,Z3,£3,5) : = 



( X, + if,i{X,,Z,,s,,d) ^ 
r3 

Z3 + S3tf/2(.X3,Zi,S3,6) 

S3 



(3.75) 
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such that Eq. (l3.70l) is transformed into 



X3 = -3X3 + a/ - Ci((5)r3 + r^h^siX^, rj, Z3, £3, 6), 
1 ^ 

f3 = 

[3 

Z3 = -2Z3 - J -63 + Si,ri,h36(X3,r3,Z-i,S3,6), 
5 

where ipiJ^s^h^^ ~ 0^(1) and i/fi ~ Op(2) are C°° functions. 
Proof of Lemma 3.8. When = 0, Eq. (|3.70l) is written as 



(3.76) 



^3 = -3^3 + ^1 -Z3 + h32(X3, 0, Z3, £3, 6), 



Z3 = -2Z3 - y + S3h33(X3,Z3,£3,S), 

5 

^3 = -2^3. 



(3.77) 



Since eigenvalues of the Jacobian matrix at the origin of the right hand side of the above are 
-3,-2,-5/2 and satisfy the non-resonance condition, there exists a C°° transformation of the 
form (^3, 23,^3) 1-^ (X3-l-(/ri(X3,Z3,e3,5),Z3-l-0'2(^3,-Z3,e3,5),e3) such that Eg. (13. 771) is linearized 
(see Chow, Li and Wang dH). The i/i'2 is of the form 1/^2 = £34^2, where 1/^2 is a C°° function, because 
if £3 = 0, Eg. (13.771) gives ^3 = -2z3 and it follows that Z3 = Z3 when £3 = 0. This transformation 
brings Eg. dXTOl) into Eg. (l3J6l) . ■ 

Proof of Prop.3.7. Note that even in the new coordinates (X3, r3, Z3, £3), the sections S3" and S""' 
are included in the hyperplanes {£3 = P3} and {r3 = pi}, respectively. 

Let us calculate the transition time T from S3" to S3"'. Since r3{t) = r3{0)e'^^ and £3(0 = 
£3(0)e^^'^^ from Eg. (|3.76l) . T is given by 



r = log 



Pi 

^3(0) 



(3.78) 



By integrating the third eguation of Eg. (13. 761) . Z3{t) is calculated as 

Z3(0 = Z3(0)e-2' + yf6p3(e-'"^ - e~^') + e''-' f p3r3(0)/z36(X3(5), 6)ds, (3. 

Jo 



79) 



where X3{s) = (,X3{s), r3{s), Z3{s), £3(5)). Owing to the mean value theorem, there exists < t = 
T(t) < t such that Eg. ( 13. 791 ) is rewritten as 



Z3(0 = (Z3(0) - V6p3)e'" + <ep3e~''l^ + p3r3{0)e-^'h3(,{X3{T),6)t. 



(3.80) 
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This and Eq. (l3.78l) are put together to yield 



Z,{T) = (Z3(0) - V6p3)(^j + V6p3(^j +P3^/l36(X3(T(r)),5)r. (3.81) 

Next, let us estimate Xt,{T). Since (X3, r3)-plane is invariant, the unstable manifold of the 
origin is included in this plane and given as a graph of the C°° function 



Xs = (f>(r,,6) = --ci(6)r, + 0(ri). 



(3.82) 



To measure the distance between ^3(0 and the unstable manifold, put X3 = (p{r^, 6) + u. Then, 
the first equation of (13.761) is rewritten as 

M = (-3 + h^jiu, r3, Z3, £3, 6))u + Z^h^^iu, r^, Z3, £3, 6) + e^h^<){u, r^, Z3, £3, 6), 

where hi,-] ~ Op{\) and h-i^, /139 are C°° functions. This is integrated as 



u{t) = e'^'Eit)\u{0)+ I e''Eisy\Z3(s)h38iuis),5) + S3is)h39(uis),d))ds 



Jo 



(3.83) 



where h( 5) = (u(s),r3(s),Z3(s),S3(s)) and E(t) = exp[j^h37(u(s),6)ds]. Substituting Eq. (l3.80l) 
and £3(0 = p^e'^'^^ and estimating with the aid of the mean value theorem, one can verify that 
u{T) is of the form 

u(T) = r3(0)X(X3(0),r3(0),Z3(0),p3,Pi,5), (3.84) 

where /z4o is bounded as r^iO) (the factor Z3(5) in Eq. (l3.83l) yields the factor r^iO)"^, and 
other terms are of O{rT,{0)^ log r3(0)). Since the transition time T is not C°° in pi and r3(0), /z4o is 
only in ^3(0), Z3(0),p3 and 6. Thus the transition map tl'^' from I^" to 2f ' along the flow of 
Eq. (l3.76l) is given by 



n 



loc 



^X3^ 
Z3 

P3 ) 



(p{pi,6) + r^h4o(X3,r3,Z3,p3,pu6) 



Pi 

v5 „5 



(Z3-V6p3)|^j +V6p3|^j -2p3^l0g|^j/l4l(X3,r3,Z3,P3,Pl,5) 



■ ^3 

P3 — 
Pi 



(3.85) 

where /74i(X3,r3,Z3,p3,pi,5) = h36(Xs(T(T)),6) is bounded as r3 ^ because X3(t(T)) is 
bounded. Since the transition time T is not C°° in pi and r3(0), /z4i is C"" in X3(0),Z3(0),p3 
and 6. Now Eq. (l3.74l) is verified by calculating O o n'°' o 0"^ Note that /?3 in Eq. (l3.74D is inde- 
pendent of r3 and pi because it comes from the inverse of the transformation (13.751) . which is of 
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the form 



' X3 +;e5(X3,Z3,P3,^) 

Z3 +P3/33(%,Z3,P3,5) 
P3 



with C°° functions fi-i and ySj. The unstable manifold fix{px,S) in (X3,r3, 23,^3) coordinate is 
obtained from that in (X3, r3, Z3, £3) coordinate as fi\{p\,5) = (p{p\,S) + il/\{(p{p\,5), 0, 0, 5). This 
proves Prop. 3.7 (I). To prove (II) of Prop. 3. 7, note that the hyperplane {r3 = 0} is invariant and 
included in the stable manifold of the origin. Since a point (^3, r3,z3,p3) converges to the stable 
manifold as r-i 0, li^^'^ix-i, ri,,z3,P3) converges to the unstable manifold as r3 ^ on account 
of the /l-lemma. This proves Prop. 3.7 (II). ■ 



3.6 Proof of Theorem 3.2 



We are now in a position to prove Theorem 3.2. Let t.y : (x, r, z, s) ^ {x - V2/3, r, z, s) be the 
translation in the x direction introduced in Sec. 3. 5. Eg. (13. 81) is obtained by writing out the map 
Tx o Tlj^ o o K23 o IIj"' o K12 o n'"^' and blowing it down to the {X, Y, Z) coordinates. At 



^loc 

first, H'S"^ o K12 o n''"^ is calculated as 



ttIoc 







r VI +^1 ] 


( 


yi 








Pi 


1 — > 


1/5 -1/5 


1 — > 




Pie/ P2 








P2 


V 



P2^^Vi+pr'% 

P-^V2+Pf^yi 

P2"^^ 
Pl^y^ 



3/5 1 



' Px+Hiip^ ' ipi+P^' Xi -qx,P2 V2+P2 J'l -qy,P2,PlSi ,P3,S)'' 

-2/5 
P3 

, T7 / -3/5 , -3/5 V -2/5 , -2/5 V 1/5 c^ 

P^ + H2(fi2 <Pi +P2 -qx,P2 <P2 +P2 Yi -qy,p2,PlSy ,P3,S) 

Pl^l 



(3.86) 



where (pi = (fiipis/ ,P2,5), V2 = V2(Pie/ Pj ,P2,S), andXi, Yi are defined by Eq. (l3.47D . 
In what follows, we omit the arguments of //i and H2. The last term in the above is further 
mapped to 



/ „3/5 3/5 TJ \ 

P3 Px + P3 tl\ 



3/5 J 



■l/5„l/5 

^4/5 „ . ^4/5 rj 
P3 Pz + P3 ^2 



P1P3 ' 



Pa 



3/5, 



3/5, 



-l/5„l/5 
'1 

4/5 , 4/5 TJ 
P3 Pz +P3 ^2 



P1P3 



P3 





X3 








Z3 




.P3. 



(3.87) 
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Let us denote the resultant as (^3, r3, Z3,p3) as above. Then, ft,^^^, proves to be given by 



n 



loc 



pi 



Pi 



(Z3 - V6p3 +P3;e3(X3,Z3,p3,(5))|^j + -logr^ ■fi4{X3,r^,Z3,P3,Pl,3) 



P3 



4/5 



+ O(eilogei) 



(3.88) 



' - V273 +/3i(pi,S)+ p\p^ '*^^e}^^/32(^3 , P1P3 ^^^e) ''^ Z3 , p3, p 1 , 5) ' 

Pi 

(Z3 - V6p3 +P3ye3(X3,Z3,p3,5)) — 

\P3 

Si 

By using the definition of in (13.291) . the third component of the above is calculated as 

(Z3- V6p3 +P3;S3(X3,Z3,P3,5)) — +0(eilogei) 

\P3/ 

= (^l + 0{p3) + H2{X,Y,p2,Pis'J\p3,5)+pfl53{x3,Z3,P3,5))sf + 0{s,\ogs,\ (3.89) 



where 



-0 -3/5 , -3/5 V V -2/5 , -2/5 v 

X = P2 ' (^1 + P2 ' Xi - q^, Y = p^ ' (p2+ P2 Yi - qy. 



From Eqs. (II3Il) and (IXTO . Eq. (IXH91) is rewritten as 

(ri + 0(p3) + mx, Y,p2) + 0(pf ) +pf Afe,Z3,P3,5))ef + 0{s, log£i). (3.90) 

Since n,'^^(xi,yi,pi, is independent of p3, which is introduced to define the intermediate sec- 
tions Ej"' and Z3", all terms including p3 are canceled out and Eg. (13. 901) has to be of the form 



(Q. + H2{X, F,P2))ef + 0{Sy log^i). 



(3.91) 



Now we look into X and Y. Since ip\{r\,e\,6) and ip2{r\,£\,5) give the graph of the center 
manifold W^{S) and since the orbit y of the first Painleve equation is attached on the edge of 
W^\5), xi = <p\{0, S\,S) and y\ = ^2(0, £i,S) coincide with y written in the ^1 coordinates. Thus 
we obtain 

O -3/5 ^ 1/5 -1/5 c^ , -3/5^ 

X = P2 ' (fiipis^' P2 ,P2,5)+P2 ^1 -q.x 

= (pf'cpi(0,p2,6) - q,) +pf'X, +pf"0{{e,lp2f") 
= pf'X,+pf'0{{s,lp2)'") 



p,"'D, [^fj' exp[-^] +P2^^^0((^,/P2)^^^) 



= D,ef'oxp[-^]+p,"'0{{sJp2y'^ 



(3.92) 
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The Y is calculated in the same manner. Functions D\,D2 and d are expanded as Eqs . (13 .48113 .491) . 
and H2 is expanded as (13.421) . Since Eg. (13. 911 ) should be independent of p2, which is introduced 
to define the intermediate sections Zj"' and Zj", Eg. (13. 911) is rewritten as 



(Q + 4(^,J/))£f + logs, 



). 



(3.93) 



where 



V A . -3/5 r d{pu5)6^ 1/ A / o -2/5 r d{pu5)6^ 
X = Di{xi,yuPi,euS)e^ exp[ J, = D2{xuyuP\,e\,S)ey exp[ J. 

(3.94) 

Similarly, since the first component of Eg. (13. 881) is independent of p2 and p3, we find that it is 
expressed as 



V273 +A(Pi,5) + J/,Pi,5)er + 0{s, logei) 



(3.95) 



with some C°° function G. 

Our final task is to blow down Eg. dX^ with Egs. dlMl [3.951 ) to the {X, Y,Z) coordinates to 
obtain Eg. (13. 81) . By the transformation (13.131) . a point (X, y,p|, e) in (X, Y, Z, e)-space is mapped 
to the point {Xp\^, Yp\^,pi,sp\^) in ^i-space. Further, it is mapped by the transition map fl^^^ to 

' - V273 + ySi (pi , 5) + G{X, J/, pi , 5)e4/5pi 4 + 0{e log e) \ 

Pi 

[q. + ti2{X, J/)) e4/5p-4 + 



in ^3 -space, in which 



r rf(pi,(5)d-| 

exp[ — ^J, 



X = D,{Xp-\Yp\\pusp-,\5)s~"'p\r-'^ ^^""^^^^ 

epi 

^(Pi,^)^ 
ep[5 



J/ = D2{Xp\\ Ypf,pu£pf, 5)e-'-l'p] exp[- 



1- 



Finally, it is blown down by (13.151) as 



(-yl2/3p] +I3i(pu6)p] + p-'G(X,}f,puS)sf + Oisi logs,)] 
Q + 4(^,J/))£f + 0(£i logei) 



By changing the definitions of Di, D2 and d appropriately, we obtain Theorem 3.2 (I) with 



Gi = - ^|2i3p\ +fi,ipuS)pl G2 = p\'G, H = H2 
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Theorem 3.2 (II) follows from the fact that the unstable manifold described in Prop. 3.7 (II) 
coincides with the heteroclinic orbit a'^{6) if written in the (X, Y,Z) coordinates. Theorem 3.2 
(III) follows from Lemma 3.4, and (IV) follows from Eg. (13. 501) because si in Eg. (13. 501) is now 
replaced by sp^^. This complete the proof of Theorem 3.2 ■ 

4 Global analysis and the proof of main theorems 

In this section, we construct a global Poincare map by combining a succession of transition maps 
(see Fig HI) and prove Theorems 1,2 and 3. 

4.1 Global coordinate 

Let us introduce a global coordinate to calculate the global Poincare map. In what follows, we 
suppose without loss of generality that the branch S ^ and 5 of the critical manifold are convex 
downward and upward, respectively, as is shown in Fig[TJ Recall that (X,Y,Z) coordinate is 
defined near the fold point and that the sections E^^, and S^^^^y are defined in Eg. (13.71) . We define 
a global coordinate transformation {x,y,z) ^ iX,Y,Z) satisfying following: We suppose that 
in the (X, Y,Z) coordinate, L+(5) = (0,0,0), L-{6) = (0,yo,Zo) with > 0,zo > 0, and that Y 
coordinates of 5^7 are larger than those of S'^, just as shown in FigHT] Let z\ > zo be a number 
and put Z2 = p\ + e~^^^ . Define the new section 

2:; = {z = pt + e-i/^'}, (4.1) 

which lies slightly above S^^. Change the coordinates so that the segment of in the region 
Z2 < Z < zi is expressed as 

{X = 0, y = -77, Z2<Z<zi}, (4.2) 

where 77 is a sufficiently small positive constant (if pi is sufficiently small). We can define such 
a coordinate without changing the local coordinate near and the expression of 11^^^ given in 
Eg. (13. 81) by using a partition of unity. We can change the coordinates near 5^7 U {Lr} in a similar 
manner without changing the coordinate expression near 5^ U {L^}. Let 

' X = /i(X,y,Z,e,5), 

- Y = f2{X,Y,Z,s,6), (4.3) 
^Z = sg{X,Y,Z,s,5\ 

be the system (11.81) written in the resultant coordinate, where the definitions of f\ , /2 and g are 
accordingly changed. 
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Fig. 1 1 : Coordinate for calculating the global Poincare map, and a slow manifold corre- 
sponding to the segment S^fe, Zi) = {(X, 7, Z) e 5^ | Z2 ^ Z < zi) of the critical manifold. 

4.2 Flow near the slow manifold 

Put 5+(z2,Zi) = {(X, Y,Z) e S^,\z2 < Z < zy}. Then, 5+(z2,Zi) is a compact attracting normally 
hyperbolic invariant manifold of the unperturbed system of (14.31) . see FigHT] In this subsection, 
we construct an approximate flow around the slow manifold corresponding to S'^{z2,Zi). 
If the parameter 5 is a constant, the existence of the slow manifold immediately follows from 
Fenichel's theorem: 

Theorem (Fenichel [U). 

Let Nhea C manifold (r > 1), and X'XN) the set of C vector fields on with the topol- 
ogy. Let F be a C vector field on A'^ and suppose that M c A'' is a compact normally hyperbolic 
F-invariant manifold. Then, there exists a neighborhood 1/ c X''{N) of the origin such that if s is 
a small positive number so that sG e U for a given vector field G 6 X''(N), then the vector field 
F + eG has a locally invariant manifold Me within an ^-neighborhood of M. It is diff"eomorphic 
to M and has the same stability as that of M. 

Further, Fenichel [|9l HOl proved that admits a fibration: there exists a family of smooth 
manifolds {'FE{p)}peM^ such that 

(i) if p ^ p', then Tsip) n Tsip') = 0. 

(ii) r,ip)nM, = {p}. 

(iii) the family {Tsip)} is invariant in the sense that (pti^eip)) c Te{4>t{p)), where 0, is a flow 
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generated by F + sG e X'iN). 

(iv) there exist C > 0,A > such that for q e T^ip), W^tip) - ^t(Q)\\ < Ce~'^', where we suppose 
for simplicity that M (and thus M^) is attracting. 

See also Wiggins ||35]| for Fenichel theory. These theorems are applied to fast-slow systems by 
Fenichel [jTTll to obtain a slow manifold M^- and a flow around M^. Roughly speaking, these 
theorems state that for a fast-slow system, there is a locally invariant manifold M^, called the 
slow manifold, within an e-neighborhood of the critical manifold M if e > is sufficiently small. 
A flow near is given as the sum of the slow motion (dynamics on M^) and the fast motion. 
If Mi; is attracting, the fast motion decays exponentially to zero and eventually a flow is well 
approximated by the dynamics on M^. 

Applying these results to our fast-slow system (14.31) . when 6 is independent of s, we obtain 
an attracting slow manifold and we can construct an approximate flow around M^. However, 
if 6 depends on e, Fenichel theory is no longer applicable in general even if e « 6. To see this, 
let us recall how the existence of Me is proved. 

For simplicity of exposition, suppose that vector fields are defined on R'" x R". We denote 
a point on this space as {x,z) e R" x R". Suppose that a given unperturbed vector field F has 
an attracting compact normally hyperbolic invariant manifold M on the subspace {x = 0}. We 
denote a flow (pt generated by the perturbed vector field F -i- eG by 



0,(X, Z, S) = ((pt (.X, Z, S), 0, {X, Z, £)). 



From the assumption of normal hyperbolicity, we can show that there exists a positive constant 
T such that 



dx 



(0,z,0) 



1 



-^(0,z,Or' < ^, for (0,z)eM, 
oz 4 



because d(p]ldx decays faster than dcf^Jdz- Since M c {jc = 0} is F-invariant, we have 



4(0,z,0) = 0, 



d(p\ 



(0,z,0) = 0, for(0,z)6M. 



(4.4) 



(4.5) 



Since the flow is continuous with respect to x, z and s, for given small positive numbers rji and 
772, there exist sq > and an open set V d M such that the inequalities 



dx 



ix,z, e) 



d(pl 



dz 



{x, z, sY 



1 



< 



\\^j(x,Z,s)\\ < 7]u 
1 

T 



dz 



{x,z, e) 



(4.6) 
(4.7) 
(4.8) 



hold for < e < So and (x, z) 6 V. Let S be the set of Lipschitz functions from M into the .x:-space 
with a suitable norm. Let 5 c be the subset of S consisting of functions h such that {h{z),z) e V 
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and their Lipschitz constants are smaller than some constant C > 0. We now define the map 
G : Sc ^ S through 



(Gh){4{h{z),z,s)) = (t>\{h{z\z,s). 



By using inequalities (I4.6[ 14.71 14.81) (and several inequalities which trivially follow from com- 
pactness of M), we can show that G is a contraction map from S c into S c- See Lemma 3.2.9 
of Wiggins [[35l . in which all inequalities for proving Fenichel's theorem are collected. Thus G 
has a fixed point hs satisfying h£(^^(hE(z),z, e)) = (l)j{hi;{z),z, s). This proves that the graph of 
X = hsiz), which defines Ms, is invariant under the flow 4>t{- , ■ , s)- The existence of a fibration 
{!^e(p)}peMs can be proved in a similar manner. 

If the unperturbed vector field F = smoothly depends on 6 and if 6 depends on s, the 
above discussion is not valid even if s « 5. The inequality (14.41) for Fg does not imply the 
inequality (14.61) for Fg + sG in general. For example, consider the linear system jc = AqX + 5AiX 
with matrices 



Suppose that 6 = ^fs. Eigenvalues of Aq+SAi are given by -S (double root), so that the derivative 
of the flow at the origin is exponentially small for t > 0. Next, add the perturbation sAjX to this 
system, where 




Although sAj is quite smaller than Aq + 6Ai if e is sufficiently small, the eigenvalues of Aq + 
6Ai + eA^ are 5 and -35, so that the derivative of the flow of the perturbed system diverges as 
t ^ oo. This shows that Eq. (l4.4l) does not imply Eq. (l4.6l) in general if 5 depends on e. Further, 
the open set Y above also depends on e through 5 and it may shrink as e ^ 0. For this linear 
system, it is easy to see that such a stability change does not occur if Aq has no Jordan block. For 
our fast-slow system, the assumption (C5) allows us to prove that such a stability change does 
not occur. 

Lemma 4.1. Let A(5, z) and 5(e, 5,z) be 2x2 matrices which are C°° in their arguments. Suppose 
that eigenvalues of A{6, z) are given by -d/j.(z, S) ± V--T<^(z, 6) with the conditions yu(z, 6) > and 
Ci>(z, 6) for 6 >0. Further suppose that 6 depends on e as e ~ o(6) (that is, s « S as s ^ 0). 
Then, eigenvalues of A(6, z) + sB(s, 6, z) are given by 




6/x(z, 6) + yTlcoiz, 6) + 0(e) 



(4.9) 



as e ^ 0. 



Proof. Straightforward calculation. 
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Now we return to our fast-slow system (14.31) . Put X = {X, Y), f = (/i, /2) and rewrite Eg. (14. 31) 



as 



X = f(X,Z,s,6), Z = sg(X,Z,s,6). 
The flow generated by this system is denoted as 



(4.10) 



4>,{X, Z, s, 6) = Z, £, 6), <Pt(X, Z, s, 6)). 



(4.11) 



Recall that 5^^(z2, Zi) is expressed as X = 0, 7 = -77; that is, /(O, -77, Z, 0,6) = for zi < Z < zi- 
When s = 0, 02(X,Z,O,5) = Z, which proves that \\{d(l)^(X,Z,0,6)/dZr^\\ = 1. Next, the 
derivative of 0j satisfies the variational equation 

-^(X,Z,0,6) = -^(0 (X,Z,0,5),Z,0,<5)^(X,Z,0,5). 
at oX oX oX 

On 5^^(z2,Zi), this is reduced to the autonomous system 

:r ^(0, -77, Z, 0, 6) = -^(0, -77, Z, 0, ^)^(0, -77, Z, 0, 6). 
dt oX oX oX 

The assumption (C5) implies that the eigenvalues of the matrix ^i^(0' -ri,Z,0,6) are given by 
-6iu\z,6) + yrTaj^(z,6). Thus 



d(f>l „ 
^{Q,-n,Z,Q,6)-0{e-'') 



on 5^(z2,2i)- This proves the inequality 



dX 



(0, -77, Z, 0,5) 



d(pl 



dZ 



(0, -77, Z, 0,5)- 



< 



4' 



(4.12) 



for some large T > 0. In general, this does not imply Eq. (l4.6l) as was explained. However, in our 
situation, by applying Lemma 4.1 to 

A(5,z) = ^(0,-77, Z, 0,5), 

it turns out that eigenvalues of the matrix t— (0, -77, Z, s, S) are of the form (14.91) for small s > 

0. Therefore, — — (0, -77,Z, e, 5) also decays with the rate 0(e ) on S^(z2,Zi). Further, the 

assumption (C6) proves that there exists a neighborhood V* of S^(z2,Zi), which is independent 
of 6, such that real parts of eigenvalues of df/dX are also of order 0{-6) on This yields 
the inequality (|4.6I) on Inequalities (|4.7|) and (|4.8I) are easily obtained. In this manner, all 
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inequalities for proving Fenichel's theorem are obtained, and the existence of the slow manifold 
and a fibration on for our system are proved in the standard way as long as s « 6 (To 

prove Theorem 3, we will suppose that 6 ~ 0(e(- loge)'''^) » s). Note that the existence of 

a neighborhood of the critical manifold, on which eigenvalues of df/dX have negative real 

parts, are also assumed in the classical approach for singular perturbed problems to estimate the 

dynamics of fast motion, see O'Malley Il271l and Smith OTI . 

We have seen that a solution of (14.101) on is written as the sum of the slow motion on the 

slow manifold and the fast motion which decays exponentially. To calculate them, it is convenient 

to introduce the slow time scale by r = st, which provides 

dX ^ dZ 
s—=fiX,Z,s,6), — = g(X,Z,s,6). (4.13) 
dr dr 

A solution of this system is given by 

X(t, €, S) = X,(t, £,6)+X f(T, £, S), 

Z(t, e, 6) = z,(r, s,6) + z /(r, s,6), 

where Xs,Zs describe the slow motion and Xf,Zf describe the fast motion. They are C°° in s 
(see Fenichel llTTl ) and their expansions with respect to s are obtained step by step according to 
O'Malley lETll as follows: We expand them as 



k=0 k=0 

OO CO 

ZsiT, e:,6) = Y, ^'zf\T, 6), Zfir, e, 5) = s'zfir, 6), 



k=0 k=0 

with the initial condition 

Z(0, s, 6) = XoiS) + 0(e), Z(0, s, 6) = zoiS) + 0(e), 

m v. At first, and zf^ are determined to satisfy the system (14.131) for e = 0. Thus is 
given by jcf ^ = (0, -t]) and zf^ is given as the solution of the equation 

dz^"^ 

-^=g(0,-r],z'^\0,S) (4.15) 
dr 

with the initial condition z[^\0, 6) = Zq(S). This system is called the slow system. Next, from the 
system (14.101) for e = 0, we obtain z^"^ = 0, and x''^ is governed by the system 

dxf dX dx^^^ 

= —{t, 0, 6) - -^(T, 6) = /((O, -77) + xf, (r), 0, 5) (4.16) 



44 



with the initial condition 

jcf (0, 6) = Xo(6) - jc,(0, 0, 6) = Xo{6) - (0, -77). (4.17) 

Fenichel's theorem (Part (iv) above) shows that if JC^^O, 6) e V^, then x''^ decays exponentially 
as ? ^ 00. In the classical approach (23, the existence of V'^ is used to estimate Eq. (l4.16l) 
directly to prove that x^j ^ decays exponentially, see also Smith OTl . To investigate behavior of a 
solution as e ^ 0, we rewrite Eq. (l4.16l) as 

-^(0,-;7,4V),0,5)xf + -ei(xf ,5), (4.18) 



where Qi ~ 0((xyf) is a C~ function. 



X 



(0)x2 

/ ^ 

Lemma 4.2. A solution of the system (14.181) is given by 

' K^{T, s) cos[i W(t)] + K2ir, e) sin[i W^(t)] K,{t, s) cos[i W^(t)] + K,(t, s) sin[i W(t)] 
^ KsiT, s) cos[iW^(T)] + Ke{T, s) sm[^W{T)] K,{t, e) co^[\W{t)\ + K^{t, e) sin[i W(t)] ^ 

exp[-- J ti\zf{s), 6)ds\{xf{0, 6) + h(t, e, 6; xf(0, 6))), (4. 19) 

where W{t) = £to^{z?\s),6)ds, Ki{i = I,-- - ,8) are C~ functions, and u ~ O{xf(0,6f) 
denotes higher order terms with respect to the initial value. 

Proof. We use the WKB analysis. Put x^p = (vi, V2) and 



^(o,-,.4"'(r),o,*) = i (4.20) 



Let us consider the linearized system 



^(^.).|<o,-,..»,o,«(:;).(:««;))(:;). 



Then, Vi(r) proves to satisfy the equation 



b'\ I ab' 
e\': -\E{a + d) + E^—\v\+\ad-bc + e{— - a) 1 vi = 0. (4.22) 



bj ' \ ' b 

We construct a formal solution of the form 



1 °° 

i(r) = exp[-_^£«S„(T)]. 



45 



Substituting it into Eq. (l4.22l) . we obtain an equation of 5o 

(S'of - (a + d)So + (ad - be) = 0. 



This is solved as 



JA+(s)ds, I A^(s)ds, 
Jo 



where 



are eigenvalues of the matrix (14.201) . For each £A+{s)ds and £A^{s)ds, 5 1, ^2, • • • are uniquely 
determined. Thus a general solution vi(t) is given by 



Vl 



1 1 

(r) = Ci exp[- J A+(s)ds^Kn(T, s) + C2 exp[- J A^{s)ds^Ku{T, s). 



where Ci,C2 e C and ^11,^12 are C°° functions. In a similar manner, it turns out that V2 is 
expressed as 

1 1 r"^ 

V2(t) = Ci exp[- J A+{s)ds\K2\{T, s) + C2 exp[- J /l_(5)J5]i^22(T, e). 

Therefore, a general solution of the system (14.211) is written as 

'exp[i j^^A^{s)ds\Kn{T,s) exp[^ J^'"/l_(^)j5]ifi2(T, e) ' 
^ exp[i j^A^{s)ds\K2i{T, s) exp[^ J^'"/l^(5)(i5]is:22(T, e) 

The fundamental matrix of (14.211) is given by 

' exp[i ^ A^{s)ds\Kn{T, s) exp[i ^ A^{s)ds\Kn{T, s) 
^ exp[^ ^ A^{s)ds\K2i{T, s) exp[^ A-{s)ds\K22{T, s) 

This shows that each component of the fundamental matrix is a linear combination of 

exp[— I //+(zf (5),(5)J5]cos[-W^(t)] and exp[— I (5), 5)^5] sin[-W^(T)]. 

^ ^0 ^ ^ Jo ^ 

Finally, the variation-of-constants formula is applied to the nonlinear system (14.181) to prove 
Lemma 4.2. ■ 



Ci 
C2 



.K2ii0,s) K22i0,s), 
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With this jc^^-*, the zeroth order approximate solution is constructed as 



Z{t, e, 6) 



zf\r,6) + 0{s),) 



0{s) 
-7] + 0(e) 



^ xfiT,6) + 0(s) 
0{s) 



(4.23) 



as long as the orbit is in V^. The first term in the right hand side denotes the position on 
and the second term denotes the deviation from Mg.. It is known that all terms x''p,z'p in the 
expansions of the fast motion decay exponentially as well as x^p ( IfTTl 17711311 ). 

Combining this approximate solution near the slow manifold with the transition map near the 
fold point, Theorem 1 is easily proved. 

Proof of Theorem 1. To prove Theorem 1 , (5 is assumed to be fixed. For the system (12.11) . take an 
initial value in V^. Then, a solution is given by (14.231) with (14.191) . These expressions show that 
when t > 0, the solution lies sufficiently close to the critical manifold 5^ if s is sufficiently small. 
Because of the assumption (A3), Zs decreases (where we suppose that S"^ is convex downward) 
with the velocity of order s (with respect to the original time scale t). Thus the solution reaches 
the section Z^, after some time, which is of order 0(1 /s). The intersection point is mapped into 
'^out by the transition map U^^^. given in Thm.3.2, and it proves that after passing through the 
distance between the solution and the orbit is of order 0{s'^^^). ■ 



4.3 Global Poincare map 

In Sec. 3, the transition map 11^^^^ around the fold point L'^(6) had been constructed. The transition 
map around the fold point L'i6) is obtained in the same way. The sections and are defined 
in a similar way to 1,^^ and (see FigjH), respectively, and the transition map 11,^^^. from an open 
set in Z7 into S^Hf along the flow of (14.31) proves to take the same form as n^^_.., although functions 
Gi, G2 and higher order terms denoted as O(eloge) may be different from one another (note that 
Q. and H are common for 11^^ and 11^"^ because they arise from the first Painleve equation). 

Since the unperturbed system has a heteroclinic orbit a' connecting L (5) with a point on 
S'^i(6) and since S'^{6) has an attraction basin which is independent of 6, there is an open set 
t/,7„y c which is independent of 6 and s, such that orbits of (14.31) starting from U'^^f go into 
V'^ and are eventually approximated by Eg. (14. 231) . Let zo be the Z coordinate of L~(6). Define 
the section to be 

Z;, = Vn {(X, Y,Z)\Y = -77, |Z - zol < Pa), (4.24) 

where p4 is a small positive number so that a solution of (14.31) starting from intersects SJ^ 
only once (see FigfT2l). 

The global Poincare map is constructed as follows: Let n^^oMr' ^hr i b^ transition maps 
from U~^^^ c into into 2^, into Z^^, respectively. Then, the transition map 11"^ from 

Kut into is given by 
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Fig. 12: The sections Z|, Sj!^ and an orbit of Eq. (l4.3l) . 



The transition map 11" from an open set in 2^^^, into is calculated in a similar manner and it has 
the same form as 11"^. The global Poincare map is given by 11+ o n . However, it is sufficient to 
investigate one of them by identifying and 2^^,^. If n+ : U~^,, 2^^,,^ is a contraction map, so 
is H"^ o n~, and if n"^ has a horseshoe, so is H"^ o H" because H"^ and H" have the same properties. 
To identify two sections and S^^,, recall that L~ = (0, 3^0. Zo) in the (X, Y, Z)-coordinate, and 
define E^,,, to be {F = - Pj}. Let i7~,, be an open set in S",^ such that the transition map 
^Hout ■ ^out ~^ ^7/ well-defined. The set i7^„, includes the point Z^^^ n a'. We identify U'^^ 
with an open set in 1,^^, by the translation 





(X] 




' X ' 


T : 




1— > 


yo -Pi 




UJ 







(4.25) 



Then, the transition map ftj^ ^^^^^ from U^^^^f c into is obtained by combining the translation 
and n|^o„f Since the velocity in the Z direction is of order s, it is expressed as 





(X\ 




{x\ 






^^11, out 


p\ 




p\ 




-77 












,Z + zo + 0(e), 



(4.26) 



where P+ is a C"" function. Since nf/o„, is C°°, we expand it as 



n 



lI,out 



(x\ 




' p{S) + o{x,z,sy 


p\ 






UJ 




, Z + zo + Ois) , 



(4.27) 
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To prove Theorem 3, we will use the fact that there exists a positive constant po > such that 
\p(6)\ > Po for < 6 < 6o, which is proved as follows: Since 6 controls the strength of the 
stability of S^, if 6 is sufficiently small, orbits which converge to (0, -?7,zo) (the intersection 
of the heteroclinic orbit or" and 5^) rotate around this point so many times. In particular, they 
intersect with SJ^ before reaching (0, -?7, Zo)- If p(d) were zero, the right hand side above tends 
to (0, -rj, zo) as X,Z,s ^ 0, which yields a contradiction. 

Next thing to do is to combine the above n//o„; with n|^^. By Eq. (l4.23l) . the transition map 
nj"„ from into SJ^ is given by 





' X ' 




' 0{s) ' 










-77 + 0(e) 






[2] 




^2 





xf(T(X,Z,e,S),S) + 0{e) 




(4.28) 



where = jc^°^(t, S) is given by (14.191) with the initial condition JC^°^(0, S) = (X, 0), Z2 = p'l + 

e'^^"^ is the Z coordinate of the section as defined before, and t = t(X, Z, e, S) is a transition 
time (with respect to the slow time scale) from a point (X,-ri,Z) to This transition time 
T is determined as follows: Let z*"\t, 5) be a solution of Eg. (14. 151) with the initial condition 
zf\0, S) = Z. Then, Eg. (14. 231) implies that r = t{X, Z, s, 6) is given as a root of the eguation 

Z2 = zf\T,6) + 0{s). 

Let t be a root of the eguation Z2 = zf\T, 6). By virtue of the implicit function theorem, r is 
written as t = t + 0(e). Since Eg. (14. 151) is independent of X and e, so is f . Thus we obtain 



T(X,Z,e,S) = f{Z,d) + 0{s). 



(4.29) 



Further, f is bounded as 6 ^ because ^ on 5^^ uniformly in < 6 < Sq. Therefore, 11^ ^ 
proves to be of the form 





' X ' 




' 0(s) " 








-77 + 0(e) 




[z 1 




Z2 



X (Kyif, e) cos[i W(f)] + K2{f, e) sin[iTO]) exp[-f £p^{zTis), 6)ds]{l + 0(e,X)) ^ 
+ X (K,{f, e) cos[i W(f )] + Kdf, e) sin[i W(t)]) exp[-f £p^(^\s), 6)ds](l + 0{s,X)) 



(4.30) 

The first line denotes the intersection point M^. n 1,'^ and thus it is independent of X and Z. The 
second line denotes the deviation from the intersection. Note that the transition map nt^ ^ from 
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into is 0{e ''''^^)-close to the identity map. Thus ^ o o nj^^^^j o 7" is calculated as 













-77 + 0(e) 


UJ 




p1 



+ 



m,I 1,11 1 1, out 



'p{6) (K,(f, s) cos[i W(t)] + K2(f, e) sin[iiy(f )]) exp[-f £iu^izf\s), 6)ds]il + 0(s, X, Z)) ^ 

p(6) (K,(f, s) cos[i W(t)] + K,(f, s) sin[iTO]) exp[-f £fi^(zT(s), 6)ds](l + 0(s, X, Z)) 



(4.31) 

where f = t(Z + zo, <5) and \t) is a solution of (14.151) satisfying the initial condition zf\0) = 
Z + zo- Finally, the transition map 

loc in,I 1,11 lI,out ' 

from f/~j^j into is obtained by combining the above map with 11^"^^. 

At this stage, we can prove Theorem 2. 
Proof of Theorem 2. To prove Theorem 2, it is sufficient to show that the map n+ has a hyper- 
bolically stable fixed point. Then, the global Poincare map (without identifying and 2^„,) has 
the same property because TV takes the same form as 11^. Indeed, if s is sufficiently small for 
fixed 5, Them. 3. 2 and Eg. (14. 311) show that the image of the map 11^ is exponentially small, and 
thus is a contraction map. Further, eigenvalues of the derivative of is of order 0{e'^^^), 
which proves that H"^ has a hyperbolically stable fixed point. ■ 



4.4 Derivative of the transition map 

If 6 is fixed, it is obvious that the transition map 11^ is of order 0{e~^^^) as £ ^ 0. However, when 
6 is small as well as s, the action of 11+ becomes more complex. In what follows, we suppose 
that 6 depends on s and s ~ o{6) (s « 5) as e ^ 0. A straightforward calculation shows that the 
derivative of 11+ is of the form 



5n+ 



diX,Z) 



X 



where (z = 1 , ■ • ■ , 4) are bounded as e — 
L5 ~o(l)asX,Z,e^ 0. 

Eigenvalues of the derivative are given by 



Li(X, Z, £, S)s^'^ L2(X, Z, £, 6)s-^'^ 
UiX, Z, £, UiX, Z, e, 6)s-^'^ 

exp[-J(p, 6)^] ■ exp[-^ ^ fi^z^Ks), 6)ds](\ + Ls(X, Z, e, 6)), (4.32) 

0, and L5 denotes higher order terms such that 



^1 = L4£-4/5 exp[-d(p, 6)^] • exp[-^ H^(zf\s), 6)ds](l + o(l)), (4.33) 
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and 



L1L4 - L2L3 1 ,5 
: e 



cxp[-d{p, 6)-] ■ exp[-- (5), S)ds](l + o(l)). (4.34) 



If 6 is fixed, they are exponentially small as e ^ 0, although if 6 is small as well as s, \Ai \ may 
become large. For example, if6 = Cs{- log sY^^ with a positive constant C, and if L4(X, Z, e, 6) 
0, is of order g-4/5g-c(-ioge)'/-^ which is larger than 1 if e is sufficiently small. On the other 
hand, I/I2I is always smaller than 1. The function L4 is given by 

OA OA oZ 

|-i^i(f, s) sin[iw^(f )] + Kjif, s) cos[ V(t)] j , (4.35) 

in which arguments of D, = D,( ■ , ■ ,pi,s,6) are given by the first and second components of 
Eq. dCTT) . From Thm.3.2 (III) and (IV), we obtain dH/dX + 0, dD^ldX ^ 0. The value p{S) is 
also not zero as was explained above. Recall that t(Z + zo, S) is defined as a transition time along 
the flow of Eg. (14. 151) . Since g < uniformly on 5^^ and < 5 < ^o, t is monotonically increas- 
ing with respect to Z. Further, iy(f ) is monotonically decreasing or monotonically increasing 
because <jj" i^^ uniformly. This proves dW(j)ldZ 4^ 0. Therefore, L4 = if and only if 

-Kx{t(Z + zo, S), s) sin[i W(t(Z + zo, 5))] + K^ir^Z + zo, S), e) cos[ V(f (Z + zo, 5))] 
= -i^i(t(zo, 5\ e) sin[^ W(t(Z + Zo, S))\ + i^2(t(zo, 5\ e) cos[i W(f (Z + zo, S))\ + 0(Z) 

is zero. If there exists Z such that the above value is zero, then it is zero for a countable set of 
values of Z because of the periodicity. For these "bad" Z, Ax degenerates and \Xx \ may become 
smaller than 1 . Now we have the same situation as the proof of the existence of chaos in Sil- 
nikov's systems. In the proof of Silnikov's chaos, an eigenvalue of a transition map degenerates 
if and only if an expression k\ sin(log(z/e)) + k2 cos(log(z/e)) is zero, where k\ and ^2 are some 
constants, see Wiggins [|34ll . 



4.5 Proof of Theorem 3 

Now we are in a position to prove Theorem 3. The proof is done in the same way as the proof 
of Silnikov's chaos. At first, we show that the transition map 11+ has a topological horseshoe: 
We show that an image of a rectangle under 11+ becomes a ring-shaped area and it appropriately 
intersects with the rectangle. Next, to prove that the horseshoe is hyperbolic, we investigate the 
derivative of n+. We can avoid "bad" Z, at which the derivative degenerates, because they are at 
most countable. 

Proof of Thin.3. Suppose that 5 = C\s{- log e)^^^ with some positive constant C\. Recall that 



51 



there exists a slow manifold within an s neighborhood of 5^^. Since it is one dimension, the slow 
manifold is a solution orbit of the system (14.101) . By virtue of Thm.3.2, this orbit intersects with 
near . Let Q e be the intersection point of this orbit and Take a rectangle R on 
Sqk? including the point Q, whose boundaries are parallel to the X axis and the Z axis (see FigjH). 
Let = C2S be the height of R, where C2 is a positive constant to be determined. The image of 
R under the map n^„(// = n^uj // ^ deformed rectangle whose "height" is also of order 0(s) 
since dZ/dt ~ 0(e). 

Next thing to consider is the shape of Il^j j o n^„f //(^). It is easy to show by using Eg. (14. 301) 
that the image of IV^^^jj{R) under the map Il^j j becomes a ring-shaped area whose radius is of 
order e'^^'^. Since the "height" of n^„,//(^) is of order s, the rotation angle of the ring-shaped 
area is estimated as 

-W{fiZ + zo + 0(s))) - -W(f(Z + zo)) = - oj\z[^\s),6)ds - 0(1). (4.36) 

^ e e Jf(z+zo) 

Thus we can choose C2 so that the rotation angle of the ring-shaped area is sufficiently close to 
2n as is shown in FigfT3l 




Fig. 13: Images of the rectangle R under a succession of transition maps. 

Finally, we consider the shape of 11^ (i?) by using Thm.3.2. Since dH/dX{0, 0) 0, the 
expansion of H is estimated as 

H{X, J/) ~ Xs~^'^ exp[-J^/e](l + 0{s^'^)). (4.37) 

This and Eg. (13. 81) show that the radius of n'^(i?) is of order Ois^^^e'^^"^). Since we put 6 = 
Cis(- log e)'''^, the ineguality 

hR = C2S « 0{s^'^e-^l') (4.38) 

holds if s is sufficiently small. Further, the ring 11^(7?) surrounds the point Q because the image 
of the rectangle R under the flow rotates around the slow manifold when passing between the 
section and St . This means that two horizontal boundaries of R intersect with the ring 11^(7?) 
as is shown in Figl5](b). It is obvious that the vertical boundaries of R are mapped to the inner and 
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outer boundaries of the ring, and the horizontal boundaries are mapped to the other boundaries 
in radial direction. This proves that the map 11^ creates a horseshoe and thus has an invariant 
Cantor set. 

To prove that this invariant set is hyperbolic, it is sufficient to show that there exist two disjoint 
rectangles Hi and H2 in R, whose horizontal boundaries are parallel to the X axis and vertical 
boundaries are included in those of R, such that the inequalities 

I|£>XII < 1' (4.39) 

immyw < i, (4.40) 

1 - ||(An2"r'|| ■ iiDXIl >2^||D,n|||-||D,n+||-||(D,np-i|p, (4.4i) 
1 - (||D,n|ll + ll(An2")-^ll) + liz).n||| • \\{D,u^,r'\\ > \\d,u^,\\ • |iAn||| • ||(D,n+r>||,(4.42) 

hold on Hi U Hj, where n| and 11:^ denote the X and Z components of 11^, respectively, and 
and D, denote the derivatives with respect to X and Z, respectively. See Wiggins [[34]| for the 
proof. We can take such Hi and H2 so that "bad" Z, at which L4 = 0, are not included. Then, 
inequalities above immediately follows from Eq. (l4.32l) : ||D;tn||| and 11^^11211 are sufficiently 
small, and ||Dj.n||| and ||Dj.n2ll are sufficiently large as e ^ 0. This proves Theorem 3. ■ 



5 Concluding remarks 

Our assumption of Bogdanov-Takens type fold points is not generic in the sense that the Jacobian 
matrix has two zero eigenvalues. However, this assumption is not essential for existence of 
periodic orbits or chaotic invariant sets. 

At first, we remark that Theorems 2 and 3 hold even if we add a small perturbation to Eq. dl.SI) . 
since hyperbolic invariant sets remain to exist under small perturbations. 

Second, we can consider the case that one of the connected components of critical manifolds 
consists of stable nodes, stable focuses and a saddle-node type fold point (i.e. a saddle-node 
bifurcation point of a unperturbed system), as in FigHH In this case. Theorem 2 is proved 
in a similar way and Theorem 3 still holds if the length of the subset of the critical manifold 
consisting of stable focuses is of order 0(1). However, analysis of saddle-node type fold points 
is well performed in |[20l l25l [T2]| and thus we do not deal with such a situation in this paper. 

We can also consider the case that one of the connected components S of critical manifolds 
has no fold points but consists of saddles with heteroclinic orbits a-, see Fig {151 In this case, 
analysis around the S is done by using the exchange lemma (see Jones IfTSlI ) and we can prove 
theorems similar to Theorems 2 and 3. Such a situation arises in an extended prey-predator 
system. In ||23]| . a periodic orbit and chaos in an extended prey-predator system are numerically 
investigated with the aid of the theory of the present paper. 
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Fig. 14: Critical manifold consisting of a saddle-node type fold point, stable nodes, and stable 
focuses and an orbit near it. 
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